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Abstract 


Near-field observations may provide tight constraints - i.e. “boundary conditions” - 
on any model of structure formation in the Universe. Detailed observational data 
have long been available for the Milky Way (e.g., Freeman & Bland-Hawthorn 2002) 
and have provided tight constraints on several Galaxy formation models (e.g.: Abadi 


et al. 2003 Bekki fe Chiba||2001 ). An implicit assumption still remains unanswered 


though: is the Milky Way a “normal” spiral? Searching for directions, it feels 
natural to look at our neighbour: Andromeda. An intriguing piece of the “puzzle” 
is provided by contrasting its stellar halo with that of our Galaxy, even more so 
since Mouhcine et al. (2005) have suggested that a correlation between stellar halo 


metallicity and galactic luminosity is in place and would leave the Milky Way halo as 
an outlier with respect to other spirals of comparable luminosities. Further questions 
hence arise: is there any stellar halo-galaxy formation symbiosis? Our first step has 
been to contrast the chemical evolution of the Milky Way with that of Andromeda by 
means of a semi-analytic model. We have then pursued a complementary approach 
through the analysis of several semi-cosmological late-type galaxy simulations which 
sample a wide variety of merging histories. We have focused on the stellar halo 
properties in the simulations at redshift zero and shown that - at any given galaxy 
luminosity - the metallicities of the stellar halos in the simulations span a range in 
excess of ~ 1 dex, a result which is strengthened by the robustness tests we have 
performed. We suggest that the underlying driver of the halo metallicity dispersion 
can be traced to the diversity of galactic mass assembly histories inherent within 
the hierarchical clustering paradigm. 
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RGB tracks (i.e. 0.8 M 0 VandenBerg et al. 2000 RGB tracks and the 
empirical NGC 6791 RGB track) which are used, as in Mouhcine et al. (2005), 
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simulation (Figures 3.2 and 3.7) at M tot 10 12 M 0 . The observa¬ 
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simulations in Renda et al. (2005b). Black dots for the gaseous par¬ 
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shown in each panel. 11601 

E.2 (continued) Galaxy rotation curves for the M to t = 5xlO n M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the 
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cosmological simulations in Renda et al. (2005b). Black dots for the 
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E.7 (continued) Galaxy rotation curves for the M tot = 10 12 M 0 semi- 


cosmological simulations in Renda et al. (2005b). Black dots for the 
gaseous particles. Red dots for the stellar particles; red boxes for 
the ensemble of stellar particles (enclosed within the magenta frame) 
which are counter-rotating with vg < 0 within each topographical 
halo at a projected radius R > 15 kpc. The size of such ensemble is 


also shown in each panel. 11661 

E.8 (continued) Galaxy rotation curves for the M to t = 10 12 M 0 semi- 


cosmological simulations in Renda et al. (2005b). Black dots for the 
gaseous particles. Red dots for the stellar particles; red boxes for 
the ensemble of stellar particles (enclosed within the magenta frame) 
which are counter-rotating with vg < 0 within each topographical 
halo at a projected radius R > 15 kpc. The size of such ensemble is 


also shown in each panel. 11671 

E.9 (continued) Galaxy rotation curves for the M tot = 10 12 M 0 semi- 
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Chapter 1 
Introduction 


“Ne osservava le rughe, 
e provava ad immaginame la storia. ” 

- Anonimo 


Near-field observations - e.g. in the Local Group and in the nearby low redshift 
Universe - of properties such as metallicity, luminosity and kinematics may provide 
tight constraints - i.e. “boundary conditions” - on any model of structure formation 
in the Universe. This “near-field” approach to cosmology is the overarching theme 
of this PhD Thesis. 

The focus narrows here on stellar halos in late-type galaxies. Understanding the 
formation history of stellar halos is one of the challenges of galactic astrophysics. 
The problem is generally framed within the context of two competing scenarios: 
one of “rapid collapse” (Eggen, Lynden-Bell & Sandage, 1962), in which the stellar 
halo is formed by the rapid collapse of a proto-galaxy within a dynamical timescale 


(~ 10 8 yr), and one of “galactic assembly” (Searle & Zinn 1978), whereby the 
stellar halo is assembled over a longer timescale (~ 10 9 yr) by the accretion of 
“building-blocks”, each with separate enrichment history. Both scenarios have their 
strengths and weaknesses, and it would appear that a hybrid model is the most 
plausible picture consistent with extant data (e.g.: Chiba & Beers 2000 Freeman & 


Bland-Hawthorn 2002). 


Detailed observational data have long been available for the Milky Way (e.g., Free- 
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man & Bland-Hawthorn 2002) and have provided tight constraints on several Galaxy 


formation models (e.g.: Abadi et al. 2003 Bekki & Chiba 2001). An implicit as¬ 
sumption still remains unanswered though: is the Milky Way a “normal” spiral? 
Searching for directions, it feels natural to look at our neighbour: Andromeda is the 
closest spiral to the Milky Way, ~ 780 kpc (more than 2 million light-years) away 
from our Galaxy. There has recently been an astonishing progress in the observa¬ 
tions of Andromeda (e.g., Ibata et al. 2007). An intriguing piece of the “puzzle” is 
provided by contrasting the stellar halo of our Galaxy with that of its neighbour, 
even more so since the Hubble Space Telescope imaging of nearby spiral galaxy 
stellar halos in Mouhcine et al. (2005) suggests that a correlation between stellar 
halo metallicity and galactic luminosity is in place and would leave the Milky Way 
halo as an outlier with respect to other spirals of comparable luminosities. Further 
questions hence arise: is there any stellar halo-galaxy formation symbiosis? Is the 
Milky Way an outlier?... Which is “normal”? 

The layout of the PhD Thesis we present here reflects the path we have under¬ 
taken to face the challenge these questions pose. Our first step has been to contrast 
the chemical evolution of the Milky Way with that of Andromeda by means of a 
semi-analytical model (Chapter [2] Renda et al. 2005"), which, in passing, has also 
led to a new understanding of the origin of Fluorine in the Milky Way (Appendix |A| 


Renda et al. 2004). We have then pursued a complementary approach through the 


analysis of several semi-cosmological late-type galaxy simulations which sample a 
wide variety of merging histories. We have first focused on the stellar halo proper¬ 
ties in the simulations at redshift zero (Chapter [3j which synthesises the background 
which is presented in Appendices |B| [H| Renda et al. 2005b), we have then traced 
galaxy mass assembly and metal enrichment in the simulations over the second half 
of the age of the Universe (Chapters [3] and [5j respectively). We finally summarise 
what (we think that) we see from our point of view in Chapter [6] 











































Chapter 2 

Contrasting the Milky Way and 
Andromeda Stellar Halos 


The chemical evolution history of a galaxy hides clues about how it formed and has 
been changing through time. We have studied the chemical evolution history of the 
Milky Way (MW) and Andromeda (M31) to find which are common features in the 
chemical evolution of disc galaxies as well as which are galaxy-dependent. We use a 
semi-analytic multi-zone chemical evolution model. Such models have succeeded in 
explaining the mean trends of the observed chemical properties in these two Local 
Group spiral galaxies with similar mass and morphology. Our results suggest that 
while the evolution of the MW and M31 shares general similarities, differences in the 
formation history are required to explain the observations in detail. In particular, 
we found that the observed higher metallicity in the M31 halo can be explained by 
either a) a higher halo star formation efficiency or b) a larger reservoir of infalling 
halo gas with a longer halo formation phase. These two different pictures would lead 
to a) a higher [O/Fe] at low metallicities or b) younger stellar populations in the 
M31 halo, respectively. Both pictures result in a more massive stellar halo in M31, 
which suggests a possible correlation between the halo metallicity and its stellar 
mass. 


2.1 Introduction 


The chemical properties of galaxies hide clues about their formation and evolution. 
Semi-analytic chemical evolution models (Talbot fe Arnett|1971 Tinsley|1980) have 
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succeeded in explaining the mean trends of galactic systems by numerically solving a 
set of equations governing the simplified evolution of the chemical elements as they 
cycle through gas and stars. One strength of these models is that they typically 
have the fewest number of free parameters, making convergence to a smaller set of 
solutions more likely. 


Strong constraints can be placed on chemical evolution models only by con¬ 
trasting them against a comprehensive set of observed properties. Since the most 


detailed observational data are generally available for the Milky Way (e.g., Freeman 


& Bland-Hawthorn 2002), successful agreement between model predictions and these 


observed properties has been obtained by several studies in the past, which help us 
to understand the formation history of the Milky Way (MW). Some models have fo¬ 
cused on the evolution of the chemical abundances both in the solar neighbourhood 
and in the whole disc, adopting a framework in which the MW has been built-up 


inside-out by means of a single accretion event (e.g.: Matteucci & Frangois 1989 


Pagel & Tautvaisiene 1995). Others have used an early infall of gas to explain the 


thick disc formation, followed by a slower infall to form the thin disc (e.g., Chi- 
appini, Matteucci & Gratton 1997). Several studies have paid particular attention 
to the chemical evolution of a larger range of elements (e.g.: Timmes, Woosley & 
Weaver 1995; Goswarni & Prantzos 2000 Alibes, Labay & Canal 2001) or focused 
on the features of the Metallicity Distribution Function (MDF) in the solar neigh¬ 
bourhood (e.g., Fenner & Gibson 2003). However, an implicit assumption remains 
unanswered: is the MW a typical spiral? 


Andromeda (M31) is the closest spiral to our MW (e.g., van den Bergh 2003). 


Previous theoretical studies of the chemical evolution of M31 have been done by Diaz 


& Tosi (1984) and Molla, Ferrini & Diaz (1996), both emphasising the evolution of 


the M31 disc. They concluded that M31 has a formation history and chemical evolu¬ 
tion similar to that of the MW. More recently, however, there has been considerable 
observational progress in the study of both galaxies. A striking difference between 
M31 and the MW is that the metallicity of the M31 halo (([Fe/H]) ~ — 0.5) []is 

significantly higher than its MW analogue (([Fe/H]) « — 1.8, e.g., Ryan & Norris 


1991), as revealed by many recent studies (e.g.: Holland, Fahlman & Richer 1996; 


Durrell, Harris, & Pritchet 2001; Sarajedini & van Duyne 2001 Ferguson & Johm 


son 


2001 Reitzel fc Guhathakurta 2002 Ferguson et al. 2002 |Worthey et al. 2005 


1 Hereafter [X/Fe] = log 10 (X/Fe) — log lo (X/Fe) 0 . 
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Bcllazzini et al. 2003 Brown et al. 2003; Reitzcl et al. 2004; Rich et al. 2004). 


It is therefore timely to attempt the construction of a chemical evolution model 
for both the MW and M31, using the same framework. Such an attempt may 
be helpful in highlighting the common features in the chemical evolution of spiral 
galaxies (at least in these two spirals) and those which remain galaxy-dependent. 


In Section 2.2 we describe our multi-zone chemical evolution model, and in 


Section 2.3 we present the results of our MW model. In Section 2.4 we show the 


results for M31. Finally, our results are discussed in Section 2.5 


2.2 The model 


In this study we use a semi-analytic multi-zone chemical evolution model for a spiral 
galaxy. This model is based on GEtool ( Fenner fe Gibson] 2003 Gibson et al.||2003 ). 
We follow the chemical evolution of both halo and disc of M31 and the MW under 
the assumption that both galaxies have formed through two phases of gas infall. The 
first infall episode corresponds to the halo build-up, and the second to the inside- 
out formation of the disc. Similar formalisms have been successful in modeling 
the chemical evolution in the solar neighbourhood (e.g.: Chiappini, Matteucci & 


Gratton 1997; Oiang et al. 1999 Alibes, Labay & Canal 2001). We assume that 
halo stars were born in a burst induced either by the collapse of a single proto- 
galactic cloud (Eggen, Lynden-Bell & Sandage 1962) or by the multiple merger of 
building-blocks (Searle & Zinn 1978 Bekki & Chiba [2001 ; Brook et al. 2003). The 
disc is assumed to be built-up inside-out by the smooth accretion of gas on a longer 
timescale. Observations of HI High Velocity Clouds (HVCs) which appear to be 
currently falling onto the MW (e.g., Putman et al. 2003 and references therein) may 
provide evidence for such gas infall. Recently, HVCs have also been detected in 
the M31 neighbourhood, though their interpretation as infalling clouds is debated 


(Thilker et al. 2004). 


We monitor the face-on projected properties of the halo and disc components. 
While this geometrical simplification is suitable for approximating the flat disc, it is 
less appropriate for the halo, whose shape is roughly spherical rather than disc-like. 
However, we consider this choice acceptable in a simplified model of the chemical 
evolution of a spiral galaxy. We follow the chemical evolution of several independent 
rings, 2 kpc wide, out to a galactocentric radius R — 10 Rj, where Rd is the disc 
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scale-length. We also ignore the bulge component, because we are interested in the 
relatively outer region (R > 4 kpc). Each ring is a single zone onto which gas falls, 
without exchange of matter between the rings. We trace the chemical evolution of 
each zone. In our model, we assume that the age of the galaxy is t now = 13 Gyr. 


The basic equations (e.g., Tinsley 1980), in a zone at a radius R, for the evolution 
of the gas surface density t) of an element i are written as follows: 


x 


x 


E gji (R,t) = 

f‘ Mb min 

ip(R,t)X i (R,t) + / i/;(R,t - T m ) 

** Mmin 

r ' M B max 


Yi(m, Z t _ Tra )—dm + k 
m 


B 


'M b 


M, 


B 


r-0.5 


f(n)ip(R,t - T m 2 )Yi(M B ,Z t _ Tm 2 )d/X dM B 




+ (1 - At) 


^(R,t - T m )Yi(m, Z t _ Trn ) 


' Me 




X 


ip(m) r ivl max 

- -dm + / ip(R,t - Tm)Yi{m, Z t _ Tm ) 

m J M Bmax 


dm + X ihalo (t)h(R)e^ 

~( f ~ t delay ) 



<p(m) 

X 


m 

+ 



^ g,i(R,t) 


s g(R,t ) 


( 2 . 1 ) 


formation rate (SFR); <p(m) is the Initial Mass Function (IMF) with mass range 
M rnm - M max ; r m is the lifetime of a star with mass m; Yj(m, Z t _ Tm ) is the stellar 


yield of the element from a star of mass m, age r m and metallicity Z t _ 


The first 


term describes the depletion of the element i which is locked-up in newly formed 
stars. The second and the fourth terms show the contribution of mass loss from 
low and intermediate mass stars. The third term describes the contribution from 
Type la SNe (SNe la). The contribution from SNe la is calculated as suggested in 
Greggio & Renzini (1983), where k, M Rrntn , M Bmax , / i min , n, /(//), r m2 are defined. 


The fifth term shows the contribution from Type II SNe (SNe II). The sixth and 
the seventh terms represent the infalling halo and disc gas, respectively. 


The Kroupa, Tout & Gilmore (1993) IMF is used here. We have chosen a lower 
mass limit of M min = 0.08 M 0 and imposed an upper mass limit of M max = 60 M Q 
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Table 2.1: The parameters for the MW and the M31 models. 



a 

^t,h (-^O j tnow) 

3d. 

kpc 

^t,d(R © Jnow) 

Th 

Gyr 

tdelay 

a d 

b d 

Vh 

v& 

M 0 pc - 2 

M 0 pc - 2 

Gyr 

MW 

2 

6 

3.0 

48 

0.1 

1.0 

2.0 

1.25 

0.125 

0.03 

M31a 

2 

6 

5.5 

46 

0.1 

1.0 

2.0 

1.25 

0.125 

0.03 

M31b 

2 

6 

5.5 

46 

0.1 

1.0 

0.7 

0.50 

12.5 

0.08 

M31c 

2 

57 

5.5 

46 

0.5 

6.0 

0.1 

0.10 

0.125 

0.08 


in order to avoid the overproduction of oxygen and recover the observed trend of 
[O/Fe] at low metallicity in the solar neighbourhood in the MW. Such IMF upper 
mass limit is currently loosely constrained by stellar formation and evolution models. 
Yet, these stellar models, and the yields they provide, are one of the most important 
features in galactic chemical evolution, although questions remain concerning the 
precise composition of stellar ejecta, due to the uncertain role played by processes 
including mass loss, rotation, fail-back, and the location of the mass cut, which 
separates the remnant from the ejected material in SNe. The SNe 11 yields are from 


Woosley & Weaver (1995). We have halved the iron yields shown in Woosley & 


Weaver (1995), as suggested by Timmes et al. (1995). The SNe la yields are from 


Iwamoto et al. (1999). The metallicity-dependent yields of Renzini & Voli (1981) 
have been used for stars in the mass range 1-8 M 0 . The lifetimes of stars as a 


function of mass and metallicity have been taken from Schaller et al. (1992). 


2.2.1 Infall 


The infall rate during the halo and disc phases is simply assumed to decline ex¬ 


ponentially, as seen by the adopted sixth and seventh terms in equation 2.1 The 
evolution of the totaj^] surface mass density, T, t (R,t ), is described by: 


/? — t t-delay) 

1 = h(R)e^ + d(R)e ~^«—, 


( 2 . 2 ) 


2 Hereafter, by “total” density we mean the sum of the stellar density and of the gas density. 
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Here, the first term describes the infall rate in the halo phase. The infall time- 
scale in the halo phase, r h , is assumed to be independent of radius, for simplicity. 
The infall of disc gas starts with a delay of t de i ay , as seen in the second term. The 
time-scale of the disc infall depends on radius as follows: 


r d (R) = a d + b a 


R 

kpc 


Gyr. 


(2.3) 


The values for the constants a d and b d are free parameters. The infall coefficients 
h(R ) and cl(R ) are chosen in order to reproduce the present-day total surface den¬ 
sity in the disc, Yt td (RR now ), and in the halo, E tt h{R,t n ow), respectively, as follows 
(e.g., Timmcs et al. 1995[): 


h(R) Yit] 1 ( y R^t now ') 


X< T h 


exp 


Th 


-1 


d(R) Yt d ( y R,t now ') 


x{^)[l-exp(- i =V_^) 


(2.4) 


(2.5) 


The infalling halo gas has been assumed to be of primordial composition. On the 
other hand, it is unlikely that the accreting gas has primordial abundance at a later 
epoch, since even low density inter-galactic medium, such as the Lyman a forest, 


has a significant amount of metals (e.g., Cowie & Songalia 1998), and it is known 
that the HVCs in the MW, which could be infalling gas clouds, have metallicities 


between 0.1 and 0.3 (Sembach et al.l 2002). Therefore, we assume that the 


gas accreting onto the disc is pre-enriched. The level of pre-enrichment can be 
loosely constrained from the observed metallicity of Galactic HVCs. We simply 
assume that the metallicity of the infalling disc material is Zi n f a u(R,t ) = Z(R,t ) 
if Z(R,t) < Z in f a u )rnax = 0.3 Zq, where Z(R,t) is the metallicity of the gas at 
the radius R and the time t, otherwise Z in f a u(R, t) = Z in f a u max = 0.3 Z Q . The 
abundance pattern of the infalling disc gas is further unknown parameter. Following 
the above simple assumption, we set the infalling disc gas, at a given galactocentric 
radius R and time t, to have the same abundance pattern as the ISM at R and t. 
This guarantees the smooth evolution of the gas abundance and of the abundance 
patterns in each radial bin. 
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2.2.2 Disc and halo surface density profiles 

We adopt the following exponential profile for the present-day total surface density 
of the disc component: 

E t , d (R,tnow ) = E ty d(R Q ,t now )e~ (R ~ Re)/Rd . (2.6) 


Here R & = 8 kpc, which is the galactocentric distance of the Sun within the MW. 
The same definition of R Q = 8 kpc is applied to the M31 models. 

The surface density profiles of the MW and M31 disc are different. We have 
chosen a scale-length of R d = 3.0 kpc for the MW (e.g.: Robin, Creze & Mohan 
1992; Ruphy et al.||1996 Freudenreich||1998 ) and R d = 5.5 kpc for M31 (Walterbos 


& Kennicutt, 1988). For the MW, we assume E t ^(.R 0 , tnow ) = 48 M 0 pc^ 2 (Kuijken 
& Gilmore 1991). We adopt E t,d(RQ,t n0 w) = 46 M 0 pc~ 2 in M31, such that the 
mass of the M31 disc is similar to that of the MW disc 10 11 M 0 as in Freeman 


1999). 


We adopt a modified Hubble law for the present-day total surface density profile 


of the halo (Binney & Tremaine 1987): 


t now ) 


J t,h 0 


1 + (. R/Rq) a 

Here we set a = 2. This corresponds to a volume halo density profile: 


(2.7) 


Pt,h{j*i tnow) 


Pt,h 0 


1 + <*> 


3/2’ 


( 2 . 8 ) 


with p tM - 2 3/ ' 1 ptj,(R Q ,t npw ) and >l, hl! - 2R Q p tM . For the MW, we assume 
Et : h(R@, t now ) = 6 M 0 pc~ 2 which yields a present-day stellar surface density at R & 
of (R e ,t now ) = 1.3 M 0 pc -2 for our model (see Section 2.3), which is consistent 
with the observed halo stellar density at the solar radius p*,/i(R 0 , tno«i) = 5.7 x 
10 4 M 0 kpc -3 as estimated by Preston, Shectman & Beers (1991). 

The assumption of a = 2, which implies oc r -3 , agrees with a recent 

analysis (Zibetti, White & Brinkmann 2003) of the halo emission for a sample of 
~ 1000 edge-on disc galaxies within the Sloan Digital Sky Survey (SDSS). This 
result is similar to the density profile oc r^ 3 ' 5 which has been suggested for the MW 


stellar halo (Chiba & Beers 2000; Chiba & Beers 2001 Sakamoto, Chiba & Beers 
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2003). 


2.2.3 Star Formation Rate 


We assume that the halo star formation (SF) happens on a short time-scale because 
of a rapid infall event associated with the collapse of a single massive proto-galactic 


cloud (Eggen et al., 1962) or with multiple mergers of building-blocks (Searle & 


Zinn 1978). The disc SF is assumed to be a more quiescent phenomenon, and likely 
to be driven by the spiral arms (e.g., Wyse & Silk 1989). Therefore, we adopt a 


different SF law for each component. 

The adopted halo SFR is described as: 


MR,*) = M ® °y r 1 p c 2 . 


(2.9) 


where zy,, is the star formation efficiency (SFE) in the halo. Therefore, the halo SFR 


follows a Schmidt law with exponent 1.5 (e.g., Kennicutt 1998). The adopted halo 


SFE is Uh = 0.125, which is approximatively half of the value = 0.25 ± 0.07) 


suggested by Kennicutt (1998) and is chosen to give the best fit to the observed halo 


MDF in the solar neighbourhood (see Section 2.3). Stars born before tdeiay , when 


the disc phase starts (Section 2.2.1), are hereafter labelled as “halo stars’ 1 


The adopted disc SFR is written as: 


Rq 


-i 


i>d(R,t) = u d E g (R,ty — M 0 Gyr * pc 


-2 


( 2 . 10 ) 


where v d is the SFE in the disc. This formulation (Wyse & Silk, 1989) reflects the 


assumption that SF in the disc is triggered by the compression of the ISM by spiral 
arms. The efficiency factor v d is a free parameter. We have found that v d affects 
both the present-day gas fraction and the disc MDF. The value v d = 0.03 is used 
in our MW model to reproduce the observed gas density profile of the MW disc and 
the observed MDF in the solar neighbourhood. 


2.3 The Milky Way model 


Using the multi-zone model described in the previous section, we construct a model 
which closely reproduces the known observational properties of the MW. The adopted 
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parameters are summarised in Table 2.1 (MW model). 


Fig. 2.1 “a” compares the model MDF in the solar neighbourhood with the 


observed MDF of K dwarfs (Kotoneva et ah, 2002). Here we have assumed a K 
dwarf mass range of 0.8 - 1.4 M 0 and convolved the model MDF with a Gaussian 
error function with a = 0.15 dex, consistent with the known empirical uncertainties 
of the observational data. There is good agreement between the model and the data. 


In Fig. 2.1 “b”, the predicted halo MDF at the solar radius is shown. Onr MDF 


includes all stars still on the main-sequence at the present-time. The MDF has been 


convolved with a Gaussian error function with a = 0.25 dex (e.g., Asplund 2003) 


Our model is in agreement with the observed halo MDF in the solar neighbourhood 


(Ryan & Norris, 1991). The evolution of [O/Fe] as a function of [Fe/H] for the gas 
component in the solar neighbourhood is shown in Fig. |2.1| “c”. The model result is 


consistent with the trend observed in local stars. Fig. 2T] “d” displays the predicted 
SF history in the solar neighbourhood. The SF history has been estimated by several 


authors (e.g.: Bertelli & Nasi 2001 Hernandez et al. 2000). Unfortunately, such 


observational estimates are not defined well enough to provide useful constraints 

“d” demonstrates that our 


on a chemical evolution model. Nevertheless, Fig. 2.1 


MW model is consistent with the broad range of the estimated SFR as summarised 
in Rana (1991). 


The left panels of Fig. 232 show the predicted radial distribution of the present- 
day hydrogen column density (upper) and of the present-day oxygen (middle) and 
sulfur (bottom) abundances of the gas, and compare them to the observations. The 
observed hydrogen surface density is obtained by summing the surface densities of 
H 2 and HI (Dame 1993) 3 The result of our MW model is compatible with the 
observed hydrogen distribution at the inner radii within the uncertainties (~ 50%) 
of the observed values Q However, the hydrogen surface density at the outer radii 
is overestimated when compared with the observations. 

The predicted present-day radial abundance profiles of oxygen and sulfur also 
reproduce the abundances observed in HII regions (Vilchez & Esteban 1996; Affler- 


bach et al. 1997), although at the outer radii the model overestimates the sulfur 


abundances. 


Our semi-analytic model, whose parameter values are summarised in Table 2.1 


3 The HI surface density profile has been recently confirmed by 


Nakanishi & Sofue 


( 2003 ). 


4 The contribution of HII should also be considered. However, so far no reliable measurement 
has been achieved. 
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Figure 2.1: Comparison between the results of the MW model (solid lines) and 
the observations of: “a” the MDF in the MW in the solar neighbourhood (closed 
boxes with error-bars, Kotoneva et ah 2002); “b” the halo MDF at R = R Q 
(Ryan & Norris 1991, closed boxes with statistical Poissonian error-bars); “c” [O/Fe] 
and [Fe/H] for stars observed in the solar neighbourhood (Carretta et al. 2000, who 
included reanalysis of Sneden et ah 1991, Tomkin et ah 1992, Kraft et ah 1992, 
Edvardsson et ah 1993; Gratton et ah 2003; Bensby, Feltzing & Lundstrom 2003; 
Cayrel et ah 2003); “d” the present-day SFR in the solar neighbourhood as sum¬ 
marised in Rana (1991, the shaded region). 
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Figure 2.2: Radial distributions of the present-day hydrogen surface density (up¬ 
per) and of the present-day oxygen (middle) and sulfur (bottom) abundances for 
the MW model (left) and the M31 models (right). In the left panels, the solid lines 
are for the MW model results. In the right panels, solid, dotted and dashed lines 
display the M3 la, M31b and M31c models, respectively. The observed distribu¬ 
tions of the total surface density of hydrogen for the MW are from Dame (1993, 
open boxes), and those for M31 are from Dame et al. (1993, open boxes) and 
Loinard et al. (1999, closed boxes). Observations of abundances in HII regions 
in the MW, from Vilchez & Esteban (1996) and Afflerbach et al. (1997) (closed and 
open boxes, respectively), are shown in the left panels. Those in M31, from Den- 
nefeld & Kunth (1981) and Blair et al. (1982) (closed and open boxes, respectively), 
are shown in the right panels. 
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satisfies the general MW observational constraints. We now use the same framework 
to study the chemical evolution of M31. 


2.4 The M31 models 


The M31 models employ for the disc a different surface density profile from that of 
the MW model. Since the halo density profile of M31 is unknown, we adopt the 


same halo profile used in the MW model (see Section 2.3). The SF history of halo 
and disc is described by the parameters listed in Table 2.1 In Section 2.4.1, we 


first show a M31 model which adopts the same parameter set of our MW model. 
This model fails to reproduce some crucial features observed in M31. Therefore, 


in Sections 2.4.2 and 2.4.3, we present two models able to better explain the key 
observations. 


2.4.1 M31a model: MW analogue 


First, we construct a M31 model (M31a) with the same parameter values of our 
MW model. These parameters are summarised in Table |2.1 The right panels of 
Fig. 2.2 show the radial distributions of the present-day hydrogen surface density 
and the radial profiles of the present-day oxygen and sulfur abundances of the gas 
phase. The radial profiles of oxygen and sulfur abundances are reproduced within the 
observational errors. The M31a model results in too high hydrogen surface density 
when compared with the data. This is because the observed hydrogen surface density 
in M31 is smaller than in the MW. 


Fig. 2.3 shows the radial profile for the mean [Fe/H] of main-sequence (MS) 
stars for the M31a model. Hereafter we simply call ([Fe/H]) the ’’mean metallicity”. 
Model results are compared with observations by Bcllazzini et ah (2003)[/ We have 
chosen as reference fields those which are mostly disc- (or halo-) dominated, with 
estimated halo- (or disc-) contamination around or lower than 10% (see Table 1 
in Bcllazzini et al. 2003): the disc-dominated fields (G287, G119, G33, G76, G322 
and G272) lie at deprojectec^] galactocentric distances R ~ 8 kpc, R ~ 12 kpc, 


5 The observed MDFs are derived from red giant branch (RGB) stars, whereas the model 
produces the MDFs of MS stars, since it is difficult to construct MDFs of RGB stars within our 
framework. We assume that this inconsistency does not significantly affect our comparison. 

6 The inclination angle of the M31 disc is *m 31 ~ 12.5°. 
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i? ~ 13 kpc, R ~ 14 kpc, R & 15 kpc and R ~ 18 kpc, respectively; the halo- 
dominated fields (G319, Gil, G351, G219 and Gl) lie at projected galactocentric 
distances R ~ 16 kpc, R ~ 17 kpc, R ~ 19 kpc, R ~ 20 kpc and R ~ 34 kpc, 


respectively. Fig. T3 shows that M3 la is in broad agreement with the observed 
mean metallicity in the disc-dominated fields, though the metallicity gradient is 


slightly steeper than the observed one (Table 2.2). On the other hand, the mean 


metallicity of the halo in M31a is too low, compared to the data. 

To overcome the failure of M31a, we explored the parameter space of the M31 mod¬ 
els which could explain the observational properties of M31, and found two viable 
solutions. In the following, we present these two models. 


2.4.2 M31b model: our best model 


The model parameters of M31b are summarised in Table 2T Compared with the 
M31a model, i.e. the MW analogue, this model has a shorter time-scale for the infall 
of disc gas (i.e. smaller ad and bj) and higher disc and halo SFE. We found that 
a combination of higher disc SFE and shorter time-scale for the disc infall leads 
to better agreement with the observed radial profile of the present-day hydrogen 
surface density. In addition, a higher SFE in the halo phase leads to a more metal- 
rich halo, and we adopt a SFE 100 times higher in the halo phase than in the 
M31a model. 

In M3 lb the radial profile of the present-day oxygen and sulfur abundances 
reproduces the observational data to roughly the same degree of the M3 la model 


(Fig. 2.2). On the other hand, the higher disc SFE leads to stronger SF and therefore 


larger depletion of gas, and lower present-day hydrogen surface density than in 
M3 la. Consequently, the result of the M3 lb model is in better agreement with the 


observed radial profile of the hydrogen surface density (Fig. 2.2). The hydrogen 


surface density profiles in both the M31b model and the observations peak at a 
radius of 8 kpc with a value of ~ 4.5 M 0 pc -2 (Loinard et ah, 1999). However, 


and similarly in the MW model, the hydrogen surface density at the outer radii is 


overestimated in M31b when compared with observations from Dame et al. (1993). 


We have tried to reproduce the observed hydrogen surface density at the outer radii 
by changing various parameters. However, we are unable to find a better parameter 
set. This might suggest that, to explain the low gas surface density at the outer 
radii in M31, another mechanism which is not included in our semi-analytic model 
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10 20 30 40 
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Figure 2.3: Radial profile of the present-day mean stellar metallicity of both disc 
(upper panel) and halo (lower panel). Solid, dotted and dashed lines represent the 
results of the M31a, M31b and M31c models, respectively. The mean metallicities 
from Bellazzini et al. (2003) are also shown (closed boxes with 1 o dispersion of their 
MDF). 
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Figure 2.4: The MDFs of the M31b and M31c models (dotted and dashed lines, 
respectively). The MDFs observed from Bcllazzini et al. (2003) are also shown 
(closed boxes). 
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is required. 

The additional benefit of the shorter time-scale of the disc infall in M31b is its 
smaller SFR at the present-day. Since M31b has a SFR which peaks at an earlier 
epoch than in M31a, M31b predicts a lower present-day SFR ( « 1.5 M 0 yr' 1 ) 
in the disc (4 - 20 kpc) than in M31a ( ~ 2.4 M 0 yr' 1 ). The observational esti¬ 
mates of the present-day SFR in M31 still have a large uncertainty. For example, 
Williams (2002) estimated the SFR for the last ~ 5 Gyr roughly between « 2 
20 M 0 yr' 1 , while Williams (2003) inferred that the total SFR for the M31 
1 M 0 yr' 1 over the past 60 Myr. This value is higher than the SFR 


and 

disc has been ~ 

of ~ 0.2 M 0 yr' 1 estimated from H a and Far Infra-Red luminosities (Devereux 


et al. 1994). In addition, all these values could be affected by systematic errors 


(e.g., Bell 2003). Nevertheless, these observations suggest that M31 has a lower 
SFR than in the MW, and M31b is consistent with this trend. 


Fig. |2.3| shows improved agreement between the results of the M31b model 
and the observations of the present-day mean metallicity for both disc- and halo- 
dominated fields. In the disc, M31b has a shallower metallicity gradient than in 
M31a, and better reproduces the observed gradient (Table [2~2| ) . This is mainly due 
to the small bd we adopted (Table [2T| . 

Fig. |2.4| directly compares the MDFs of M31b with the MDF obtained in Bellazz- 


ini et al. (2003) in different fields. The model MDFs are derived from main-sequence 


stars, and convolved with a Gaussian error function with a = 0.25 dex. There is 
qualitative agreement between the MDFs of the M3 lb model and those observed, es¬ 
pecially in the disc-dominated fields. The main exception is the G1 halo-dominated 
field, whose observed MDF is more metal-poor and bimodal. The MDFs of the G351 
and G219 halo-dominated fields show narrower peaks than those of the model. Dif¬ 
ferences in the shape of the observed halo MDFs among the different fields suggest 


possible inhomogeneities in the M31 halo (see also Bellazzini et al. 2003). 

In the framework of onr semi-analytic models, M31b gives the best fit for the 


observational constraints available in M31. In Section 2J3, we discuss the formation 
process of M31 which is implied by these results. 


2.4.3 M31c model: another possible model 

The results of the M31b model suggest that a stronger halo SF can produce a metal- 
rich halo. However, applying the higher SFE is not the only way to induce stronger 
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Table 2.2: Radial gradients (dex/kpc) for the present-day mean stellar metallicity of 
the MDFs in the M31 models. The results for the MW model and for the M31 fields 
observed in Bellazzini et al. (2003) are also presented as reference. (A([Fe/ H]) / AR) 
refers to the range 4-16 kpc for the MW model and 4-20 kpc for the M31 mod¬ 
els. (A{[Fe/H})/AR) halo refers to the range 4-30 kpc in the halo for both the 
MW and the M31 models. 



(A([Fe/H])/AR) 

(A<[Fe/tf])/Afl) Mo 

MW 

-0.058 

-0.021 

M31 obs 

-0.026 

-0.008 

M31a 

-0.041 

-0.021 

M31b 

-0.021 

-0.003 

M31c 

-0.014 

-0.009 


SF. This can also be achieved by employing a larger reservoir of infalling halo gas, 
without increasing the SFE. Here we discuss the M31c model, which assumes a much 
larger present-day halo surface density and a longer (6 Gyr) halo phase (Table [24~] ). 
Again, to explain the observed hydrogen surface density, M31c uses a shorter disc 
infall time-scale, and a weaker radial dependence, which also leads to a present-day 
SFR similar to that of the M3 lb model. 


Fig.[2[2 shows the predicted present-day radial distributions of oxygen and sulfur 
abundances which are in agreement with the observational data to roughly the same 
degree of the other models. The present-day hydrogen surface density profile of the 
M31c model is consistent with Dame et al. (1993) in the inner region of the disc at 
R < 12 kpc. However, M31c has a broader peak at a larger radius than in M31b. 
Consequently, the model overestimates the surface density in the outer region at 
R > 12 kpc. Again, this result might suggest that additional physical mechanisms 
operate, to explain the observed low hydrogen surface density at the outer radii. 


Fig. |2.3| shows that M31c also reproduces the observed present-day radial distri¬ 
butions of the metallicity, similarly to the M31b model. The metallicity gradient in 
M31c is flatter than in M3 lb at the inner radii as a consequence of the weaker radial 
dependence of the disc infall time-scale (i.e. a smaller bd). In the halo fields, M31c 
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also has systematically higher metallicity than in M31a and a steeper gradient than 
in M31b. Although the gradient is in better agreement with the observations than 


in M31b (Table 2.2), M31c has lower mean metallicities than observed, especially 
at the outer radii. In Fig. 2.4 , the predicted MDFs at different radii are displayed. 
The MDFs of M31c are also in good agreement with those observed, especially in 
the disc. On the other hand, the halo MDFs show a more pronounced metal-poor 
tail than in M31b. 

Thus, the M31c model also leads to acceptable results, although M31b is in 
better agreement with the observations. The next section discusses the formation 
history of M31 implied by the M31b and M31c models. 


2.5 Discussion 

The results in the previous section have shown that the main trends of the chemical 
evolution history of both the MW and M31 can be described using the same frame¬ 
work. We have also found that to explain the observations in more detail, different 
sets of parameter values for the formation history are required for the MW and M31, 
respectively. In this section, we discuss how these different parameters values relate 
to different galaxy formation histories. 

The previous sections have shown that the main differences between the MW 
and M31 include the present-day hydrogen surface density in the disc and, more 
significantly, the present-day halo metallicity. The smaller hydrogen surface density 
in M31 can be explained by a combination of the shorter disc infall time-scale and 
the higher disc SFE. 

The most striking difference is that the mean metallicity observed in the M31 
halo is about ten times higher than that observed in the MW halo, though both 
galaxies have similar mass and morphology. 


Which is a “typical” halo? Metal-rich halos seem more common than metal- 


poor ones, as pointed out by 

Zibetti et al. 

(2004 

) who analysed a 

sample of about 

1000 edge-on disc galaxies within the SDSS. 

Harris & Harris 

(2001 

) have also shown 


that the NGC5128 halo MDF closely resembles that in the M31 halo. This similarity 
is perhaps suggestive of a common history in the halo formation of both galaxies, 
despite their different Hubble-type. 


A straightforward way to obtain a more metal-rich halo would be by requiring ho- 
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mogeneous pre-enrichment of the infalling halo gas with a metallicity Z & 0.1 Zq. 
However, this metallicity seems too high, if such pre-enrichment occurred homoge¬ 
neously in the whole universe. For example, although this value is close to the mean 
metallicity of Damped Lyman a systems (DLAs), many DLAs have much lower 
metallicity ( Pettini| 2003 Prochaska et al.||2003 ). If DLAs are more evolved systems 
than the infalling halo gas, their metallicity should be higher than that induced by 
pre-enrichment. Therefore, a pre-enrichment of Z ^0.1 may not be tenable. 

We found that the observed metal-rich halo of M31 could be explained by two 
scenarios without any pre-enrichment: a) a higher halo SFE; b) a larger reservoir 
of infalling halo gas with a longer halo phase. Scenario a) might be explained in 
a hierarchical clustering regime, where the halo is formed by accretion of building- 
blocks. Theory predicts that gas removal by supernova-driven winds should be more 
effective in lower mass systems, leading to a consequent suppression of the SFR and 


Dekel & Silk 

1986 

Efstathiou 

2000) 


based on the SDSS also suggest that the SFE decreases with decreasing stellar mass 
in low mass galaxies with a stellar mass M* < 3 x 10 10 M 0 (Kauffmann et ah, 
2003). Therefore, higher SFE in the M31 halo phase would be explained if the 


building-blocks of the M31 halo are systematically higher mass systems than those 


of the MW halo (see also Harris & Harris 2001). This result again supports the 
notion of diversity in the formation of spiral galaxies which are apparently similar 
in mass and morphology and belong to the same Local Group. Admittedly, our 
model does not adhere to the hierarchical-clustering scenario in a self-consistent 
manner, however, our chemical evolution models can be interpreted as the “mean” 
SF and chemical evolution history of the stars which end up at each radius. 

As an observational consequence of our models, we found that both scenar¬ 
ios a) and b) predict a more massive stellar halo in M31, respectively about 6 and 
9 times more massive than in the MW, whose stellar halo mass is ~ 10 9 M 0 . This 
result suggests that there might be a correlation between the halo metallicity and 
its stellar mass. Using dynamical simulations, Bekki, Harris & Harris (2003) show 
that the stellar halo comes from the outer part of the progenitor discs when the 
bulge is formed by a major merger of two spiral galaxies. Based on this, they also 
predict the correlation of the metallicity of the stellar halos and the mass of the 
bulges which were formed by major mergers, since larger bulges have the larger pro¬ 
genitors, and progenitor spiral galaxies should follow the observed mass-metallicity 
relation. Although they do not mention a correlation between the masses of the 



















22 


CHAPTER 2. CONTRASTING THE MILKY WAY AND ANDROMEDA 


bulge and halo, it would naturally be expected. Thus, a major merger scenario 
could explain our conclusion. Future observational surveys will better quantify the 
correlation between the halo metallicity and its stellar mass (e.g., 


Mouhcine et al. 


2005). 


It is possible to distinguish scenarios a) and b) through observation. First, due 
to a longer halo phase, scenario b) predicts intermediate-age population in the M31 
halo. This picture could explain the recent evidence of intermediate-age population 


by deep imaging of the M31 halo (Brown et al. 2003). 


It is worth mentioning that the metallicity gradient in the stellar halo is sensitive 
to the assumed present-day stellar halo density profile, especially in scenario b); a 


steeper density profile leads to a steeper metallicity gradient. Pritchet fe van den 


Bergh (1994) suggested that the outer halo of M31 can be modeled by a power law 
surface brightness profile of / oc R~ 4 , which is much steeper than what we assumed 
(oc R ~ 2 ). We found that such a steep profile rules out scenario b) to reproduce the 
flat metallicity gradient observed in the M31 halo. However, it is still difficult to 
accurately measure the halo surface brightness of M31 (e.g., Zibetti et al. 2004). 


Thus, more observational estimates of the M31 halo surface brightness profile would 
be an important test for this scenario. 

We also found that the higher halo SFE in scenario a) leads to about 0.2 dex 
higher [O/Fe] when [Fe/H]< -1 dex, due to intense halo SF. Although measuring 
[O/Fe] is a hard challenge for the current available instruments, this task could be 
accomplished by the next-generation large-aperture optical telescope. 


2.5.1 Prospect 

The framework we have used can explain the main trends in the chemical properties 
of both the MW and M31. However, recent observations of stellar streams both in 


the MW (e.g.: Helmi et al. 1999 Chiba & Beers 2000 Ibata et al. 2002 Brook et 


al. 2003| |Navarro, Helmi fe Freeman 2004; Majewski et al. 2004) and M31 (Ibata 
et al. 2001| McConnachie et al. 2003 Merrett et al. 2003 Zucker et al. 2004 Lewis 


et al. 2004) clearly identify inhomogeneities in the chemical and dynamical history 


of both galaxies, which could suggest that a significant fraction of the halo stars 
results from late accretion of satellite galaxies. In light of these recent observations, 
it will be important to study both the MW and M31 in more detail by employing 


a self-consistent chemo-dynamical model (e.g., Brook et al. 2004) to trace their 




























































































2.5. DISCUSSION 


23 


interrelated chemical and dynamical histories. 
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Chapter 3 

Stellar Halo 

Metallicity—Luminosity Relationship 
for Late—Type Galaxies 


The stellar halos of late-type galaxies bear important chemo-dynamical signatures 
of galaxy formation. We present here the analysis of 89 semi-cosmological late-type 
galaxy simulations, spanning ~ 4 magnitudes in total galactic luminosity. These 
simulations sample a wide variety of merging histories and show significant dis¬ 
persion in halo metallicity at a given total luminosity - more than a factor of 10 
in metallicity. Our analysis suggests that galaxies with a more extended merging 
history possess halos which have younger and more metal rich stellar populations 
than the stellar halos associated with galaxies with a more abbreviated assembly. A 
correlation between halo metallicity and its surface brightness has also been found, 
reflecting the correlation between halo metallicity and its stellar mass. Our simu¬ 
lations are compared with recent Hubble Space Telescope observations of resolved 
stellar halos in nearby spirals. 


3.1 Introduction 


Understanding the formation history of stellar halos is one of the classical pursuits 
of galactic astronomy. The problem is generally framed within the context of two 
competing scenarios: one of “rapid collapse” (Eggen, Lynden-Bcll & Sandage, 1962), 
in which the stellar halo is formed by the rapid collapse of a proto-galaxy within 
a dynamical timescale (~ 10 8 yr), and one of “galactic assembly” (Searle & Zinn, 


25 





26 CHAPTER 3. THE HALO METALLICITY-LUMINOSITY RELATIONSHIP 


1978), whereby the stellar halo is assembled on a longer timescale ( 10 9 yr) by 


the accretion of “building-blocks”, each with separate enrichment history. Both 
scenarios have their strengths and weaknesses, and it would appear that a hybrid 


model is the most plausible option consistent with extant data (e.g.: Chiba & Beers 


2000 Freeman & Bland-Hawthorn 2002). 


An intriguing piece of the halo formation “puzzle” is provided by comparing the 
stellar halo of our own Milky Way with that of its neighbour, M31. First, despite 
their comparable total galactic luminosities, the stellar halo of M31 is apparently 
much more metal-rich than that of the Milky Way (e.g.: Ryan fe Norris|1991 Mould 


fe Kristian 1986 Durrell, Harris & Pritchet 2001; jFerguson et al. 2002 Brown et 


al. 2003 Bellazzini et al. 2003 Durrell et al. 2004 Ferguson et al. 2005 Irwin et 


al. 2005). In fact, the halo of M31 bears a closer resemblance to that of NGC 5128 


(e.g., Harris & Harris 2001), despite their differing morphological classifications. 

The stellar halo-galaxy formation symbiosis has been further brought to light 
by the recent work of Mouhcine et al. (see also Tikhonov, Galazutdinova & Droz- 
dovsky 2005). The deep Hubble Space Telescope imaging of nearby spiral galaxy 
stellar halos in Mouhcine et al. (2005) suggests a correlation between stellar halo 
metallicity and total galactic luminosity. On the surface, this correlation appears 
to leave our own Milky Way halo as an outlier with respect to other spirals, with a 
stellar halo metallicity ~ 1 dex lower than spirals of comparable luminosities (such 
as M31, as alluded to earlier). 

However, it must be noted that the metallicity of the Galactic halo in the Solar 
Neighbourhood (e.g., in Ryan & Norris 1991) comes from spectroscopic metallicities 
in a kinematically-selected sample, whereas that of the M31 halo, as well as those of 


the stellar halos of nearby spiral galaxies (Mouhcine et al., 2005), have primarily been 
derived from photometric metallicities in topographically-selected samples. Since it 
is now becoming possible to obtain spectroscopic metallicities for significant samples 


of kinematically-selected giants in the M31 halo (e.g.: Ibata et al. 2004 Kalirai et 


al. 2006 Chapman et al. |2006 Ibata et al. 2007), it will be crucial to assess the 


consistency of the M31 stellar halo metallicity as derived from kinematically- and 
from topographically-selected samples, respectively. 

A grasp of the true scatter around the general trend in the halo metallicity- 
luminosity relation awaits a larger observational data set, however any theory which 
attempts to explain such a relation needs to simultaneously account for the apparent 
metallicity discrepancy between the Milky Way and the Andromeda stellar halos. 


























































































3.2. SIMULATIONS 


27 


The question arises as to what is driving the scatter of halo metallicity for galaxies 
of comparable luminosities? 

In what follows, we investigate whether differences in the stellar halo assem¬ 
bly history can explain the diversity seen in halo metallicities. Using chemo- 
dynamical simulations, we have constructed a sample of 89 late-type galaxies, span¬ 
ning ~ 4 magnitudes in luminosity, and sampling a wide range of assembly histories 
at a given luminosity (or mass). We contrast the stellar halo metallicity-galactic 
luminosity relation in the simulations with the recent empirical determination of 


Mouhcine et al. (2005). The framework in which the simulations have been con¬ 


ducted is described in Section [3T2] , while Sections |3.3| and [T4] present the results and 
the related discussion, respectively. 


3.2 Simulations 


The simulations employed here are patterned after the semi-cosmological adiabatic 


feedback model of Brook et al. (2004), and were constructed using the cherno- 


dynamical code GCD+ (Kawata & Gibson, 2003). GCD+ self-consistently treats the 


effects of gravity, gas dynamics, radiative cooling, star formation, and chemical 
enrichment (including Type II and la Supernovae and relaxing the instantaneous 
recycling approximation). In this framework, gas within the SPH smoothing kernel 
of Type II Supernova explosions is prevented from cooling for the lifetime of the 
lowest mass star which ends as a Type II Supernova (Brook et al. 2004 and references 


therein). Details of both GCD+ and the feedback model can be found in Kawata & 


Gibson (2003) and Brook et al. (2004), respectively. 


This semi-cosmological version of GCD+ is based on the seminal work of Katz & 


Gunn (1991). The initial condition for a given simulation is an isolated sphere of dark 


matter and gas. This “top-hat” overdensity has an amplitude <5; at initial redshift z ; , 
which is approximately related to the collapse redshift z c by z c = 0.36h;(l + Zj) — 1 


(e.g., Padmanabhan 1993). We set z c = 2.0, which determines hi at Zi = 40. Small- 
scale density fluctuations based on a CDM power spectrum are superimposed on the 


sphere using C0SMICS (Bertschinger 1998), and the amplitude of the fluctuations is 


parameterised by erg. These fluctuations are the seeds for local collapse and subse¬ 
quent star formation. Solid-body rotation corresponding to a spin parameter A is 
imparted to the initial sphere to incorporate the effects of longer wavelength flue- 
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tuations. For the flat CDM model described here, the relevant parameters include 
O 0 = 1, baryon fraction fib = 0.1, H 0 = 50 km s _1 Mpc -1 , spin parameter 
A = 0.06, and erg = 0.5. We employed 14147 dark matter and 14147 gas/star par¬ 
ticles. 

A series of 89 simulations - all with the same number of particles, collapse redshift 
z c , and spin parameter A - were completed spanning a factor of 50 in mass in four 
separate mass “bins”, each with a related initial comoving radius R4 0 at z ; = 40: 
M to t = 10 11 M 0 (R4 0 = 0.7 Mpc); 5xlO n M 0 (1.1 Mpc); 10 12 M 0 (1.4 Mpc); 
5xl0 12 M 0 (2.4 Mpc). 

For each total mass, we run simulations with different patterns of small-scale 
density fluctuations which lead to different hierarchical assembly histories. This was 
controlled by setting different random seeds for the Gaussian perturbation generator 

in COSMICS. 


3.3 Results 

Our grid of 89 simulations was employed to populate the redshift z = 0 stellar 


halo metallicity-luminosity plane (i.e. ([Fe/H])haio - My) shown in Fig. 3.1 (filled 
symbols). We have been conservative in our topographical definition of “stellar 
halo”, adopting a projected cut-off radius of 15 kpc to delineate the halo from 
(possible) disc and bulge contaminants. Such a choice should also allow an easier 
comparison with the observations of halo fields in nearby spiral galaxies. 

The stellar metallicity distribution (MDF) was then generated using all the stel¬ 
lar particles from the “halo region”, convolved with a Gaussian with crpe/H] = 0.15 dex, 


representing the typical observational uncertainties (e.g., Bcllazzini et al. 2003). 
Only those simulations with > 100 halo stellar particles are included in the analysis 
here. The related MDFs are displayed in Appendix |Bj 

For each MDF, the peak is displayed in Fig. |3.1[ 68% of the stellar particles 
within a given (simulated) MDFQ typically span 1.4 dex in metallicity. We con¬ 
firmed that the halo MDFs are insensitive to projection effects, even in the most 
massive simulations, suggesting that our cut-off for the “halo region” is sufficient for 
minimising disc and bulge contamination. The optical properties of our simulated 
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Figure 3.1: The stellar halo metallicity-total galactic V-band luminosity relation 
(([Fe/H])haio ~ My). A filled box corresponds to the peak of the MDF for a simu¬ 
lation. The related MDFs are displayed in Appendix [Bj Open (unlabelled) boxes 
represent the data of Mouhcine et al. (2005); labelled boxes represent additional 
data taken from the literature (see the text for details). 68% of the stars in the 
simulated (observed) MDFs are typically enclosed within ±0.7 dex (±0.35 dex) of 
the peak of the related MDF. 





130 CHAPTER, 3. THE HALO METALLICJITY-L UMINOSITY EELATIONSHIPl 


& Lanqon (2003), taking into account the age and the metallicity of each stellar 


particle. We note that the V-band luminosity shown along the abscissa of Fig. |3.1 
is the total galactic luminosity (i.e. halo + bulge + disc). 


Fig. 3.1 also shows the observed halo metallicity-luminosity values for 13 nearby 


spirals, primarily taken from Mouhcine et al. (2005), supplemented with data from 


the literature (Binney & Merriheld 1998 Elson 1997 Harris, Harris & Poole 1999; 


Brooks, Wilson & Harris 2004(0 For the observational datasets we located the MDF 
peak exactly as we did for the simulated datasets. The related MDFs are displayed 
in Figures [~B.1| and |~B.2[ 68% of the stars in the observed MDFs are typically enclosed 
within ± 0.35 dex of the MDF peak, a factor of ~ 2 narrower than the simulation^ 
What is readily apparent from Fig. 3.1| is that significant dispersion in halo 
metallicity (> 1 dex) exists at any given total galactic luminosity in our simulations. 
For a set of simulations with the same initial total mass, the only difference among 
these runs can be traced to the random pattern of the initial small-scale density 
fluctuations; this translates directly into differing hierarchical assembly histories. 
Qualitiatively, it would appear that the assembly history alone may account for the 
diversity in halo metallicity at a given galactic luminosity, and thus account for the 
apparent outliers in the observed trend. The dispersion in halo metallicity at any 


given total galactic luminosity is summarised by the histograms in Fig. T3, where 
the distribution of the MDF peaks for the semi-cosmological simulations in Rcndaet 


al. (2005b) is shown, although such histograms should be taken as plain summaries, 


since for each simulation the pattern of the initial small-scale density fluctuations 
- and the related hierarchical assembly history - has been arbitrarily chosen rather 
than drawn from a hypothetical a priori distribution of patterns. 


Fig. |3.2| shows snapshots of the gas particles at redshift z = 1.5 and the associated 
z = 0 stellar Age Distribution Function (ADF) for three of the simulations with 
M tot = 10 12 M 0 . These representative runs show that simulations with more metal- 
rich halos are assembled over a longer timescale, and thus possess a broader ADF 
(in both the halo and the associated galaxy). Such a scenario is consistent with 


the evidence, presented by Brown et al. (2003), of a substantial intermediate-age 


2 Note that the identification of the stellar halo in M33 is still debated (e.g., Tiede, Sarajedini 
& Barker 2004; see also McConna chie et aL| 2006). 


3 Our current model does not take into account any pre-enrichment scenario due to extremely 
metal-poor stars (Pop III hereafter) whose detailed physics is still much debated (e.g.: Woosley, 


Heger & Weaver 2002; |Larson 2005). An early and homogeneous pre-enrichment of the proto- 
galactic masses due to Pop III could lead to narrower MDFs at lower redshift. 
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Figure 3.2: Upper panels: Projected distribution of the gas particles at redshift 
z = 1.5 for simulations with M tot = 10 12 M 0 . Lower panels: The associated stellar 
age distributions at z = 0 for the galaxies in the upper panels. The solid (dotted) 
histogram corresponds to the stellar age distribution for the entire galaxy (stellar 
halo). The corresponding z = 0 halo metallicities are denoted in each panel. 
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Figure 3.3: Distributions of the z = 0 MDF peaks for the semi-cosmological simu¬ 
lations in Renda et ah (2005b) with M to t = : 5xlO n M 0 ; 10 12 M 0 ; 5xl0 12 M 0 . 
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Figure 3.4: The relation between the mean stellar halo age and the peak of the halo 
MDF in our simulations. The halo ADFs are displayed in Appendix |D| 



34 CHAPTER 3. THE HALO METALLICITY-L UMINOSITY RELATIONSHIP 



<age> (Gyr) 


Figure 3.5: The relation between the mean halo age and the mean age for the entire 
galaxy in our simulations. Both the halo and the galaxy ADFs for each simulation 
are displayed in Appendix |D| 
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Figure 3.6: The relation between the peak of the halo [O/Fe] distribution and the 
peak of the halo MDF in our simulations. The related MDFs are displayed in 
Appendix [Bj The halo [O/Fe] distributions are displayed in Appendix [Cj 
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Figure 3.7: Upper panels: Stellar halo metallicity-halo I-band surface brightness 
relation (([Fe/H]) halo - /xj 13,10 ). In the left panel, our fiducial simulations with 
M to t = 10 12 M 0 and 14147x2 particles are shown, while the middle panel shows 


the corresponding lower-resolution simulations with 5575x2 (crosses) and 9171x2 
(open boxes) particles. The error bars show the 68% Confidence Level around each 
MDF peak. The average number of halo particles is shown in the bottom right 
corner of each panel. Simulations with fewer than 100 stellar halo particles (at 
R > 15 kpc) are not included here. Simulations A - C of Fig. 3.2 are displayed in 
the left panel. The right panel shows the ([Fe/H])h a io - /ii iaI ° relation for all the sim¬ 
ulations. Lower panels: Halo metallicity-stellar mass relation (([Fe/H])h a io - M^ al °). 
In the left panel, our fiducial simulations with M tot = 10 12 M 0 , while the middle 
panel shows the corresponding lower-resolution simulations. The right panel shows 
the ([Fe/H])haio ~ M* al ° relation for all the simulations. 
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metal-rich population in a topographically-selected halo field of M31. Conversely, 
simulations with more metal-poor halos are assembled earlier through more of a 
monolithic process, with a consequently narrower ADF. The relationship between 


the mean stellar halo age and the peak of the halo MDF is summarised in Fig. 3.4 


We note that the mean halo age in the simulations broadly correlates with the mean 


galaxy age as shown in Fig. 3.5, thus suggesting that the assembly history of the 
halo “knows” of the formation history of the other galaxy structural components 
(e.g., bulge and disc). Both the halo and the galaxy ADFs for each simulation are 
displayed in Appendix |Dj 


Further, the most metal-poor stellar halos (which, recall, formed preferentially 
via more of a monolithic collapse) possess ct-elements to iron ratios a factor of ~ 2 
higher than the most metal-rich halos (which formed preferentially over a more 
extended period of hierarchical clustering), as shown in Fig. 3.6 Such a trend is 
expected if the different amounts of a-elcments and Iron released by Type II and la 
Supernovae over different timescales is taken into account (e.g.: Timmes, Woosley 
& Weaver 1995; Woosley et al. 2002). The halo [O/Fe] distributions are displayed 
in Appendix [Cj 


Future analysis of higher resolution simulations will clarify whether or not the 


relationships Figures [T4l |3.5| and |3.6| point to are broadened in our sample because 
of the current resolution - thus mimicking tighter correlations. 


The upper panels of Fig. |3.7 show the stellar halo metallicity-I-band surface 
brightness relation (([Fe/H])haio - /ii ial °) for simulations with M tot = 10 12 M 0 ; the sur¬ 
face brightness was measured at a projected distance of 20 kpc from the dynamical 
centre of each simulation. The halo metallicity-stellar mass relation (([Fe/H])haio ~ 
M* aI °) is shown in the lower panels of Fig. 3.7 The three galaxies presented in 
Fig. T2 are also displayed in Fig. 3.7| An immediate correlation is apparent with 
more massive halos possessing higher surface brightness and also higher metallicity. 
This is consistent with a picture in which galaxies that experienced more extended 
assembly histories have more massive stellar halos, with both higher halo metallicity 
and halo surface brightness - i.e. higher stellar halo density. 
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3.3.1 Robustness 

Resolution 

The issue of model convergence (resolution) is always a concern when interpreting 
cosmological simulations (particularly when including baryons). To test this, we 
conducted a series of simulations for M tot = 10 12 M 0 , with 5575x2 and 9171x2 
particles, to supplement the default grid (which used 14147x2 particles); the middle 
panels of Fig. |3.7| show where these lower-resolution simulations sit in the halo 
metallicity-surface brightness plane and in the halo metallicity-stellar mass plane, 
respectively. 

The consistency between higher (left panels) and lower (middle panels) resolution 
simulations is reassuring} 


Halo Semantics 


Up to now we have dealt with a topographical definition of stellar halo: we have 
labelled as “halo” the ensemble of stellar particles in a simulation - at z = 0 - at 
a projected radius R > 15 kpc. What if the “halo” is defined in a different way? 
That’s to say: given a different “halo” labelling, how robust is the dispersion in 
halo metallicity at a given total luminosity? As a robustness test, we have tried to 
include kinematical information in our halo definition - in a simple way. We have 
relabelled as “halo” the ensemble of stellar particles in a simulation - at z = 0 - 
at a projected radius R > 15 kpc and counter-rotating i.e. with circular velocity 
vg < 0. The galactic circular velocity profiles of the sample we have analysed are 
displayed in Appendix |Ej Here we focus on the stellar halcj^} Only those simulations 
with > 100 stellar particles which are counter-rotating with vg < 0 at a projected 
radius R > 15 kpc are included in the robustness test here. The related MDFs are 
displayed in Appendix [F| The MDF peaks are shown in Fig. 

The dispersion in metallicity at a given total luminosity we found among the 
topographical halos is not wiped out when we include kinematical information in 


3.8 


4 Since 

halo 


limit our 


analysis 
-2 


to simulations with > 100 halo stellar particles, the 


we 

> 28 mag arcsec - ^ region of the halo metallicity-surface brightness plane and the 


< 10 9 M 0 region of the halo metallicity-stellar mass plane are underpopulated by the 


O 

Mhaio 

lower-resolution runs (middle panels of Fig. 3.7). 

5 Clearly, the galactic circular velocity profiles show a range of disc features among the late-type 
galaxy simulations we have analysed. Such features and their relations with both the framework 
we have chosen at the current resolution and the patterns of initial density fluctuations i.e. the 
formation histories we have arbitrarily sampled will be the topic of future studies. 
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Figure 3.8: The stellar halo metallicity-total galactic V-band luminosity relation 
(haio([F e /H])kme ~ My). A hlled box corresponds to the peak of the MDF which 
refers, for a simulation, to the ensemble of counter-rotating (vg < 0) stellar particles 
within each topographical halo (at a projected radius R > 15 kpc). The related 
MDFs are displayed in Appendix [F] The observational data label-code is as in 


Fig. 3.1 
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our halo definition, as shown in Fig. |3.8| - to be compared with Fig. m 

Given a simulation, the relationship between hak3 ([Fe/H]) topos (the “topographi¬ 
cal” halo metallicity) and halo ([Fe/H]) kine (the metallicity of the ensemble of counter¬ 


rotating stellar particles within each topographical halo) is shown in Fig. 3.9 whereas 


the discrepancy in metallicity between the two “halo” labels is shown in Fig. |3.10 
as a function of the “topographical” halo metallicity. 

For each subset, the distribution of the “topographical” halo mass fraction the 
ensemble of its counter-rotating stellar particles amounts to is shown in the right 


panels of Fig. 3.11 


Both Figures 3.9 and 3.10 support what is summarised by the histograms in the 
left panels of Fig. 3.11 i.e. the discrepancy between halo ([Fe/H]) topos and halo ([Fe/H]) kme 
is distributed around ~ 0 dex within a range of « ± 0.5 dex - although the signif¬ 
icance of such a distribution is undermined by both the small size and the current 
resolution of the sample of simulations which have been analysed. A larger sample 
size and higher resolution in the simulations will allow a stronger statistical analysis 
and more complex kinematical selections in the future. 


Halo Metallicity 

Up to now, each halo MDF we have analysed and the related metallicity peak 
have referred to the ensemble of stellar particles within what has been labelled 
as “halo” in a simulation, either the ensemble of stellar particles at a projected 
radius R > 15 kpc (the “topographical” halo) or the ensemble of counter-rotating 
(vq < 0) stellar particles within the “topographical” halo. 

As a robustness test, we generate the simulated Colour-Magnitude-Diagram 
(CMD hereafter) for each “topographical” halo. Next, we approach each simulated 


CMD the same way an observed CMD is approached in Mouhcine et al. (2005), i.e., 


through the same pipeline, a metallicity-colour relationship is constructed out of 
the fiducial Red Giant Branch (RGB hereafter) tracks to derive the metallicity of 
each (generated) star in the RGB of the (simulated) CMD. 

The Padova^ theoretical isochrones (Girardi et ah, 2002) in the metallicity range 
of 0.0001 < Z < 0.030 have been used to generate the simulated CMD out 

of the ensemble of stellar particles in the “topographical” halo for each simula¬ 
tion. The theoretical isochrones are taken into account up to the RGB tip and 


J http://pleiadi.pd.astro.it 
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Figure 3.9: The relationship between halo ([Fe/Fl]) topos (the “topographical” halo 
metallicity) and halo ([Fe/H]) kine (the metallicity of the ensemble of counter-rotating 
stellar particles within each topographical halo). 
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Figure 3.10: The discrepancy between halo ([Fe/H]) topos (the “topographical” halo 
metallicity) and halo ([Fe/H]) kine (the metallicity of the ensemble of counter-rotating 
stellar particles within each topographical halo) as a function of the “topographical” 
halo metallicity. 
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Figure 3.11: Left panels: For each subset, the distribution of the discrepancy be¬ 
tween hak) ([Fe/Fl]) topos (the “topographical” halo metallicity) and halo ([Fe/H]) kine 
(the metallicity of the ensemble of counter-rotating stellar particles within each 
topographical halo). Right panels: For each subset, the distribution of the “topo¬ 
graphical” halo mass fraction the ensemble of its counter-rotating stellar particles 
amounts to. 
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Figure 3.12: The Padova theoretical isochrones (Girardi et ah, 2002) up to the 
RGB tip in the metallicity range of 0.0001 < Z < 0.030 which have been used to 
generate the simulated CMD out of the ensemble of stellar particles in the “topo¬ 
graphical” halo for each galaxy simulation. They are shown as black dots. Red dots 
display the fiducial RGB tracks (i.e. 0.8 M 0 VandenBerg et al. 2000 RGB tracks and 
the empirical NGC 6791 RGB track) which are used, as in Mouhcine et al. (2005), to 
construct a metallicity-colour relationship to derive [Fe/H] out of the (V - I) colour 
of each (generated) star in the RGB of the (simulated) CMD at -3.5 < 1 < -3.0 
(i.e. the shaded region). The histogram in the top panel displays the (V - I) colour 
distribution for the Padova theoretical isochrones at -3.5 < 1 < -3.0. 
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are shown in Fig. |3.12| as black dots. Data storage has constrained us to choose 
5% halo mass resolution i.e. for each stellar particle a Simple Stellar Population, 
whose total mass amounts to 5% of the stellar particle mass, is generated, with a 


Salpeter (Salpeter, 1955) Initial Mass Function (IMF hereafter) over a mass range of 
0.1 < m*/M Q < m top where m top = m^ B is the largest stellar mass in the theo¬ 
retical isochrones - at the given stellar particle age and metallicity - for a star behind 
the RGB tip. As a conservative choice, neither age nor metallicity extrapolation is 
performed. Given the IMF shape and both the age and the metallicity range of the 
ensemble of stellar particles the simulated CMD is generated out of, we consider the 
information-loss such IMF mass range implies not a significant one. The observa¬ 
tional errors are taken into account (e.g., Seeker & Harris 1993) and are randomly 
drawn for I < -1 from a Normal distribution with cp = cry = exp(I — 0.6) and 
ay _ i = \j + cry which implies cp ~ 0.01 at I = -4.0 and op ~ 0.2 at I = -1 
whereas (since our focus is on the upper RGB) at I > -1 we set ep = cry = 0. 
The simulated CMDs for the simulations labelled as A - C in Figures [3.2| and |3T 


are shown in Figures 3.13 - 3.15 at a resolution (0.5%) lower than that (5%) of the 
corresponding simulated CMDs because of image storage constraints. 

Once the simulated CMD is generated, both the same fiducial RGB tracks 


(i.e. 0.8 M 0 VandenBerg et al. 2000 RGB tracks and the empirical NGC 6791 RGB track) 
and the same pipeline as in Mouhcine et al. (2005) are used to construct a metallicity- 
colour relationship to derive [Fe/H] out of the (V - I) colour of each (generated) star 
in the RGB of the (simulated) CMD at -3.5 < I < -3.0. The fiducial RGB tracks 
are shown in Fig. 3.12 as red dots. The related (V - I) colour distributions for each 
simulated CMD at-3.5 < I < -3.0 are displayed in Appendix [G] whereas the re¬ 
lated MDFs derived via metallicity-colour relationship (through the same pipeline 
Mouhcine et al. 2005) are shown in Appendix [HJ The MDF peaks are shown 


as m 


in Fig. |3.16| 

The dispersion in metallicity at a given total luminosity we found among the 
topographical halos is not wiped out when the halo metallicity is derived from the 


halo simulated CMD via metallicity-colour relationship, as shown in Fig. 3.16 - to 


be compared with Fig. 3.1 


Given a simulation, the relationship between lialo ([Fe/H]) topos (the “topographi¬ 
cal” halo metallicity) and halo ([Fe/H]) CMD (the “topographical” halo metallicity as 
derived from the halo simulated CMD via metallicity-colour relationship) is shown 
in Fig. 3.17 whereas the discrepancy in metallicity between halo ([Fe/H]) topos and 




























46 CHAPTER 3. THE HALO METALLICITY-LUMINOSITY RELATIONSHIP 



V-l 


Figure 3.13: Topographical halo simulated CMD for the label “A” semi-cosmological 
simulation (Figures |3.2| and [3.7[ ) at M tot 10 12 M 0 . The observational errors are taken 
into account (e.g., Seeker & Harris 1993) and are randomly drawn for I < -1 from 
a Normal distribution with (j\ — cr v = exp(l — 0.6) and cry - i = a/c 2 + cry 
which implies <j\ ~ 0.01 at 1 = -4.0 and cri ~ 0.2 at I = -1 whereas at I > -1 
(j i = cry = 0. Red dots display the simulated CMD without taking obser¬ 

vational errors into account. Black dots display the simulated CMD when taking 
observational errors into account. 
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V-l 


Figure 3.14: Topographical halo simulated CMD for the label “B” semi-cosmological 
simulation (Figures |3.2| and [3.7[ ) at M tot 10 12 M 0 . The observational errors are taken 
into account (e.g., Seeker & Harris 1993) and are randomly drawn for I < -1 from 
a Normal distribution with ai = cry = exp(l — 0.6) and cry _ i = \Jcr\ + cry 
which implies aj m 0.01 at 1 = -4.0 and o\ ~ 0.2 at 1 = -1 whereas at 1 > -1 
(j i = ay = 0. Red dots display the simulated CMD without taking obser¬ 

vational errors into account. Black dots display the simulated CMD when taking 
observational errors into account. 
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Figure 3.15: Topographical halo simulated CMD for the label “C” semi-cosmological 
simulation (Figures |3.2| and [3.7[ ) at M tot 10 12 M 0 . The observational errors are taken 
into account (e.g., Seeker & Harris 1993) and are randomly drawn for I < -1 from 
a Normal distribution with cp = cry = exp(l — 0.6) and cry - i = a/ct 2 + <jy 
which implies <j\ ~ 0.01 at 1 = -4.0 and cp ~ 0.2 at I = -1 whereas at 1 > -1 
cri = cry = 0. Red dots display the simulated CMD without taking obser¬ 

vational errors into account. Black dots display the simulated CMD when taking 
observational errors into account. 
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Figure 3.16: The stellar halo metallicity-total galactic V-band luminosity relation 
(haio([Fe/H])CMD _ M v ). A filled box corresponds to the peak of the MDF which is 
derived out of the (V - I) colour distribution at -3.5 < 1 < -3.0 via metallicity- 
colour relationship through the same pipeline as in Mouhcine et al. (2005). The 
metallicity of each artificial star of the ensemble the MDF refers to is derived via 
metallicity-colour relationship out of the (V - I) colour of the generated star in 
the RGB of the simulated “topographical” halo CMD at -3.5 < 1 < -3.0. The 
related (V - I) colour distributions are displayed in Appendix [G] whereas the related 
MDFs derived via metallicity-colour relationship are shown in Appendix [H} The 
observational data label-code is as in Fig. 


3.1 
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Figure 3.17: The relationship between lialo ([Fe/H]) topos (the “topographical” halo 
metallicity) and halo ([Fe/H]) CMD (the “topographical” halo metallicity as derived 
from the halo simulated CMD via metallicity-colour relationship). 
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halo 


([Fe/H]) CMD is shown in Fig. 3.18 as a function of halo ([Fe/H]) topos . 


The histograms in Fig. 3.19 display the distribution of the discrepancy between 
halo ([Fe/H]) topos and halo ([Fe/H]) CMD , for each subset, and show that halo ([Fe/H]) CMD 
is generally biased towards higher metallicities (< 0.3 dex more metal-rich as shown 


in Fig. 3.19) because the fiducial RGB tracks do not extend towards colours as blue as 


those the Padova theoretical isochrones extend towards, as shown in Fig. 3.12 which 
(given that, as a conservative choice, no extrapolation has been performed beyond 
the range covered by the fiducial tracks) amounts to narrowing the metallicity range 
which is sampled at low metallicities, once the metallicity is derived out of the 
metallicity-colour relationship. 

Interestingly, since the fiducial RGB tracks have globular cluster age, the ro¬ 
bustness test above allows us to evaluate which bias could be brought in by such 
uniformly old halo age assumption when the halo metallicity is derived out of the 
metallicity-colour relationship by interpolating among fiducial RGB tracks (e.g., as 
Mouhcine et al. 2005). Such bias can be quantified by the histograms in Fig. 3.19 


m 


even if, admittedly, the construction of each simulated CMD out of the correspond¬ 
ing “topographical” halo and its subsequent analysis via fiducial RGB tracks have 
been constrained by our conservative choice not to perform any extrapolation be¬ 
yond the range covered by either the theoretical isochrones we used to generate the 
simulated CMDs or the fiducial RGB tracks we used to analyse them. The dis¬ 
crepancy between liaIo ([Fe/H]) topos (i.e. the metallicity of the “topographical” halo 
as derived by the analysis of the MDF of the stellar particles in the simulation 
which constitute the “topographical” stellar halo ensemble - the stellar particles in 
the simulation at a projected radius R > 15 kpc) and halo ([Fe/H]) CMD (i.e. the 
metallicity of the “topographical” halo as derived from the halo simulated CMD via 
metallicity-colour relationship built out of the fiducial RGB tracks with globular 
cluster age) amounts to < 0.3 dex. Future analysis of higher resolution simulations 
would further improve such comparison. 


3.4 Discussion 

We have presented here an analysis of the characteristics of late-type galaxy stel¬ 
lar halos formed within a grid of 89 simulations, with particular emphasis placed 
upon the relationship between stellar halo metallicity and the associated galactic 
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Figure 3.18: The discrepancy in metallicity between lialo ([Fe/H]) topos (the “topo¬ 
graphical” halo metallicity) and halo ([Fe/H]) CMD (the “topographical” halo metal¬ 
licity as derived from the halo simulated CMD via metallicity-colour relationship) 
as a function of the “topographical” halo metallicity. 
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Figure 3.19: For each subset, the distribution of the discrepancy between 
halo ([Fe/H]) topos (the “topographical” halo metallicity) and halo ([Fe/H]) CMD (the 
“topographical” halo metallicity as derived from the halo simulated CMD via 
metallicity-colour relationship). 
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luminosity. It helps to stress here that “Halo Semantics” is currently a contro¬ 
versial topic (e.g., Ibata et al. 2007 and references therein): we have labelled as 
“topographical halo” the ensemble of stellar particles in a simulation at a projected 
radius R > 15 kpc, at z = 0. Although future analysis of simulations at resolutions 
higher than the current would further improve the study we have undertaken, we 
have shown that - at any given total luminosity or conversely total dynamical mass - 
the stellar halo metallicities in the simulations span a range in excess of ~ 1 dex, a 
result which is strengthened by the robustness tests we have performed. 

We suggest that the underlying driver of this metallicity dispersion can be traced 
to the diversity of galactic mass assembly histories inherent within the hierarchical 
clustering paradigm. Galaxies with a more protracted assembly history possess more 
metal-rich and younger stellar halos, with an associated greater dispersion in age, 
than galaxies which experience more of a monolithic collapse. 

For a given total luminosity (or dynamical mass), those galaxies with more ex¬ 
tended assembly histories also possess more massive stellar halos, which in turn leads 
to a direct correlation between the stellar halo metallicity and its surface brightness 


(as anticipated by earlier semi-analytical models - e.g., Renda et al. 2005). By ex¬ 
tension, such a correlation may prove to be a useful diagnostic tool for disentangling 
the formation history of late-type galaxies. 


Recently, Mouhcine et al. (2005) have presented an observed correlation between 

The 


stellar halo metallicity and total galactic luminosity, as shown in Fig. 3.1 


observed dispersion in the mean halo metallicity at a given galactic luminosity is 
smaller than what we find in our simulations. However the latter can account 
for the outliers in the observed trend. Since our motivation has been to study 
which is the effect of the pattern of the initial density fluctuations alone on the 
stellar halo features at redshift z = 0 in late-type galaxy simulations, it helps to 
note that galaxy formation, as it is observed, is an ongoing process which is the 
result of the interplay among different parameters, of which the pattern of initial 
density fluctuations (thus the merging history) is one. We have shown that the 
merging history alone may be held responsible of the dispersion in halo metallicity 
at comparable total galactic luminosities, as apparently observed for example in our 


Milky Way and in Andromeda (see Section 3.1). 

This begs the question... Which is normal ? Our study suggests that if the stellar 
halo was (primarily) assembled through more of a monolithic collapse, such a low 
metallicity is indeed what should be expected; conversely, the fact that the M31 
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stellar halo is significantly metal-rich is suggestive of a more protracted assembly 
history (see also Hammer et al. 2007). An observational consequence of these dif¬ 
fering formation histories is the prediction that the M31 stellar halo should possess 
a high surface brightness; observations tentatively support this prediction (Reitzel, 
Guhathakurta & Gould 1998; Irwin et al. 2005). 

Further observations are needed to tighten our grasp of the strength and scatter 
of the stellar halo metallicity-luminosity relation, which we have shown to be a 
useful diagnostic tool for disentangling the formation history of late-type galaxies. 
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Chapter 4 

Galaxy Mass Assembly since z ~ 1 


We analyse 112 N-body/hydrodynamical late-type galaxy simulations to study the 
stellar mass assembly and the relationships between stellar, gaseous and total mass 
in galaxies since z ~ 1 down to the present day, with an emphasis on the effects of 
the merging histories on the z = 0 galaxy properties. 

The redshift evolution of the stellar mass in simulations with comparable lumi¬ 
nosities at z = 0 shows a significant dispersion, which lies in their differing merging 
histories. This contrasts with the redshift evolution of the total mass within the 
central 100 kpc, which is already assembled by z ~ 1.5. The stellar-to-total mass 
ratio in the simulations broadly agrees with the observed (e.g., Conselice et al. 2005). 
We find that massive simulations are generally more evolved than their lower mass 
counterparts. 


Conselice et al. 2005 


4.1 Introduction 


The hierarchical structure assembly has become a well established paradigm in the 
last decade. The properties of the cosmic microwave background and the large 


scale structure (e.g.: Spergel et al. 2003 Bahcall et al. 1999) support a picture 
where a large fraction of matter is non-baryonic, gravitationally interacting only, 
and assembled in a hierarchical fashion. Galaxy formation in such a scenario is a 
continuous process where galaxy properties are the results of its merging history, 
mass and environment. 

Observationally, much progress has been made in studying the global star for¬ 
mation history and the build-up of the stellar mass over a significant fraction of 
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the Hubble time (e.g.: Lilly et al. 1996 Madau et al. 1996). However, many details 
are still missing. The analysis of galaxy stellar and total masses at different cosmic 
epochs can thus assist in constraining the likely scenarios of galaxy evolution. 


Consclice et al. (2005) have recently found a lack of evolution in the galaxy 


stellar-to-total mass ratio over the redshift range 0.2 < z < 1.2 (see also Bohm & 


Ziegler 2006), thus implying that if galaxies continue to build up over the second 


half of the Hubble time stellar and dark components are then growing together. 

Here, we focus on galaxy mass assembly over the last half of the age of the 
Universe. The results of the simulations, the comparison with the data, and the 


implications of our results are described in Sect. |4.2| We summarise our conclusions 
in Sect. 14.31 


4.2 Results 


We have constructed an ensemble of 112 late-type galaxy simulations in a chemo- 
dynamical framework (Section |3.2[ ), spanning a factor of ~ 50 in total mass, and 
sampling a range of assembly histories at a given total mass. 

This sample includes the same 89 semi-cosmological GCD+ simulations analysed 
in Renda et al. (2005b) and 23 more simulations, completed within the same frame¬ 


work, with the same parameters described in Section 322, and distributed among 
the same mass “bins”. These latter simulations sample more patterns of small-scale 
density fluctuations, which lead to different hierarchical assembly histories. A sub¬ 
set of 9 simulations (with M to t = lO n M 0 , collapse redshift z c = 1.5 and spin 
parameter A = 0.054) samples more density fluctuations patterns. 

As a benchmark for the semi-cosmological framework, we have also analysed 


2 among the disc galaxy simulations presented in Bailin et al. (2005), which are 


fully cosmological GCD+ runs and adopt a A-dominated CDM cosmology (Bailin et 


al. 2005 and references therein). The photometric properties of a stellar particle in a 
simulation are modeled as those of a Simple Stellar Population (SSP hereafter) with 
the same age and metallicity. The evolution of the SSP properties as a function of 


age and metallicity is taken from Mouhcine & Langon (2003) 


In the following, the galaxy total mass, including both dark and baryonic (i.e. gaseous 
and stellar) matter, is measured for the region within 100 kpc from the centre of the 
stellar mass distribution, at each redshift, whereas the stellar and the gaseous mass 
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are measured within 15 kpc from the centre of the stellar mass, at each redshift. 


Figures 4.1, 4.2 and 4.3 display the redshift evolution of the total, stellar and 
gaseous mass normalised to the z = 0 value for each simulation. Fig. |4.4| displays 
the redshift evolution of the stellar-to-total mass ratio whereas Fig. |4.5| displays the 
redshift evolution of the gaseous-to-stellar mass ratio. Black (red) lines show the 
semi-cosmological (cosmological) runs. Magenta for the semi-cosmological simula¬ 
tions with collapse redshift z c = 1.5. 

The consistency between the assembly histories of total, stellar and gaseous mass 
in both semi-cosmological and cosmological simulations is reassuring as a robustness 
test for the semi-cosmological framework. 

The assembly of the total mass within the inner 100 kpc in the simulations is 
~ completed by z > 1.5, independently of their total mass. On the other hand, 
the stellar mass assembly histories within the central 15 kpc are more extended, 
i.e. at z ~ 1, ~ half of the baryonic mass in the inner region is gaseous, with a 
large variety of assembly patterns. Massive simulations have assembled their stellar 
content early on (typically by z ~ 1.5) and have burnt a significant fraction of their 
gas supply by now, others continue to form stars to much later epochs, whereas the 
simulations with collapse redshift z c = 1.5 have their baryonic content dominated 
by the gaseous component at all redshifts. Galaxies with comparable stellar masses 
can have significantly different assembly histories of their stellar content, and diverse 
gaseous contents. This suggests that galaxy stellar mass should be considered neither 
a robust tracer of the merging history nor a good tracer of the total baryonic galaxy 
mass. 

Although the slope of the baryonic-to-total mass ratio does not significantly 
evolve over the redshift range which is analysed (as displayed in Fig. 4.6), Fig¬ 
ures 4.7 and 4.8 show that, in our sample, galaxies of any given stellar mass, and 


at any given redshift, span a wide range of gaseous-to-stellar mass ratios. Star 
formation does not uniformly proceed in all galaxies: some galaxies have converted 
a significant fraction of their gas reservoirs into stars, whereas in others the baryonic 
content is still dominated by the gaseous component. 

The relationship between stellar mass (within the 15 kpc central region) and total 
mass (within the 100 kpc central region) in the simulations is shown in Fig. 4.9, at 
different redshifts. Stellar and total mass in the simulations are correlated since 
z ~ 1 down to the present day, thus suggesting that galaxy stellar and dark masses 
grow together continuously, similarly to the results drawn from the observations in 
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0 0.5 1 1.5 0 0.5 1 1.5 

redshift z redshift z 


Figure 4.1: The redshift evolution of the total mass within the 100 kpc central 
region, normalised to the z = 0 value for each simulation. The semi-cosmological 
simulations with collapse redshift z c = 2 are shown as black lines. Magenta for the 
collapse redshift z c = 1.5 runs. The cosmological simulations are shown in red. 
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0.5 1 1.5 
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Figure 4.2: The redshift evolution of the stellar mass within the 15 kpc central 
region, normalised to the z = 0 value for each simulation. The semi-cosmological 
simulations with collapse redshift z c = 2 are shown as black lines. Magenta for the 
collapse redshift z c = 1.5 runs. The cosmological simulations are shown in red. 
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Figure 4.3: The redshift evolution of the gaseous mass within the 15 kpc central 
region, normalised to the z = 0 value for each simulation. The semi-cosmological 
simulations with collapse redshift z c = 2 are shown as black lines. Magenta for the 
collapse redshift z c = 1.5 runs. The cosmological simulations are shown in red. 
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Figure 4.4: The redshift evolution of the stellar (within the 15 kpc central region) 
to total (within the 100 kpc central region) mass ratio, normalised to the z = 0 
value for each simulation. The semi-cosmological simulations with collapse redshift 
z c = 2 are shown as black lines. Magenta for the collapse redshift z c = 1.5 runs. 
The cosmological simulations are shown in red. 
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0 0.5 1 1.5 0 0.5 1 1.5 

redshift z redshift z 


Figure 4.5: The redshift evolution of the gaseous-to-stcllar mass ratio (within the 
15 kpc central region), normalised to the z = 0 value for each simulation. The semi- 
cosmological simulations with collapse redshift z c = 2 are shown as black lines. 
Magenta for the collapse redshift z c = 1.5 runs. The cosmological simulations are 
shown in red. 
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Figure 4.6: The relationship between baryonic mass (i.e. gaseous and stellar mass 
within the 15 kpc central region) and total mass (within the 100 kpc central re¬ 
gion), at different redshifts. The semi-cosmological simulations with collapse red- 
shift z c = 2 are shown as filled symbols (heptagons for M tot = 5 x 10 12 M 0 , 
hexagons for M tot = 10 12 M 0 , pentagons for M tot = 5 x 10 11 M 0 , boxes for 
M tot = 10 11 M 0 ). Larger empty boxes for the collapse redshift z c = 1.5 runs. 
The cosmological simulations are shown as large open triangles. 
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Figure 4.7: The relationship between gaseous mass and stellar mass (within the 
15 kpc central region), at different redshifts. Symbols are the same as in Fig. 4.6[ 
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Figure 4.8: The relationship between baryonic gas fraction and stellar mass, at 


different redshifts. Symbols are the same as in Fig. 4.6 
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Figure 4.9: Stellar mass (within the 15 kpc central region) against total mass (within 
the 100 kpc central region), at different redshifts. Symbols are the same as in 
Fig. 4.6 The stellar-to-total mass ratio in the simulations is compared with the 


ratio drawn from the observations in Conselice et al. (2005) and Bohm & Ziegler 
(2006), displayed as six and four vertices stars, respectively. 
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Conselice et al. 

(2005 

) and 

Bohm & Ziegler 

(2006) 


The scatter around the stcllar-to-total mass ratio in the simulations increases 
with redshift. However, it is smaller than the dispersion around the observed rela¬ 
tionship. The simulations analysed here are sampling neither a full range nor an 


a priori distribution of parameters (Section 3.2). Varying such parameters might 
increase the scatter in the simulations. Also, the dispersion around the observed 
relationship could be partially due to systematics in the observations. In fact, it 
is difficult to accurately estimate the dynamical mass for an intermediate redshift 
galaxy. The dynamical masses of these objects have been usually estimated using 


circular velocities obtained by long slit spectroscopy. However, Flores et al. (2006) 
have recently pointed out that long slit spectroscopy reveals overdispersed galaxy 
properties, because it tends to bias the circular velocity in galaxies with disturbed 
velocity fields. 


As previously shown in Fig. 4.6 a tighter correlation exists between total mass 
and baryonic mass. Also, the baryonic-to-total mass ratio is less dispersed than 
the stellar-to-total mass ratio, similarly to what is observed for z ~ 2 galaxies 


(Erb et ah, 2006), thus reinforcing that the galaxy stellar mass should be considered 
neither a robust tracer of the galaxy merging history nor a good tracer of the total 
baryonic galaxy mass. 


4.3 Conclusions 

We have presented the analysis of a sample of semi-cosmological late-type galaxy 
simulations with stellar mass larger than ~ 10 9 M 0 , to explore whether hierarchical 
formation scenarios may account for the constraints on galaxy mass assembly over 
the second half of the age of the Universe. 

The main results can be summarised as follows: 

- Stellar, gaseous, and total masses are found to be correlated in the simulations 
since z ~ 1 down to z = 0, in broad agreement with the scaling relations 
observed in local and intermediate redshift galaxies. 

- Galaxies with comparable stellar masses in the simulations can have signifi¬ 
cantly different assembly histories of their stellar content, and diverse gaseous 
contents. This suggests that the galaxy stellar mass should be considered nei- 
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ther a robust tracer of the galaxy merging history nor a good tracer of the 
total baryonic galaxy mass. 



Chapter 5 

The Mass—Metallicity Relation 
since z ~ 1. 


We analyse 112 N-body/hydrodynamical late-type galaxy simulations to study the 
metal enrichment of the stellar component and the interstellar medium since z ~ 1, 
down to z = 0. The relationships between stellar mass and metallicity for both 
stellar and gaseous components in the simulations at z = 0 are in broad agreement 
with the relationships locally observed. The dispersion around these relationships 
in the simulations, which lies in their diverse merging histories, broadly agrees with 
the observed. We find that the integrated stellar populations in the simulations are 
dominated by stars as old as 4 - 10 Gyr. For massive simulations, the age range 
of the integrated stellar populations agrees with the observations. On the contrary, 
simulations with stellar mass ~ 10 9 M 0 at z = 0 tend to be older than locally 
observed galaxies with comparable stellar masses. 


5.1 Introduction 

The correlation between galaxy metallicity and luminosity in the local universe is 


one of the most significant observational results in galaxy evolution. Lcqueux etal 


(1979) first revealed that Oxygen abundance increases with total mass for irregular 


galaxies. The luminosity-metallicity relation for irregulars was later confirmed by 
Skillman et al. (19891) among others. Subsequent studies have extended the relation 


to spiral galaxies (e.g.: Garnett & Shields 1987 Zaritsky et al. 1994 Garnett et al. 


1997), and to elliptical galaxies (Brodie & Huchra 1991). More recently, large sam- 
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pies of star-forming galaxies drawn from galaxy redshift surveys, e.g. the 2dF Galaxy 
Redshift Survey and the Sloan Digital Sky Survey (SDSS hereafter), have been used 


to confirm the lunrinosity-metallicity relation over a broader range (Lamareillc et 


al. 2004 Tremonti et al. 2004). The stellar populations of local galaxies with low 


stellar mass are found to be generally young and metal-poor, whereas massive galax¬ 


ies are found to be old and metal-rich (e.g., Gallazzi et ah 2005). Different groups 
have explored the properties of the interstellar gas at earlier epochs in intermediate- 


Hammer et ah 

2001 

Lilly et al. 2003 

Liang et ah 

2004 

Maier et ah 2004 


Kobulnicky & Kewlcy 2004; Maier et ah 2005 Mouhcine et ah 

2006 

Lamareillc et 

ah 

2001 

) and high-redshift galaxies (1.5 < z < 4: 

Pettini et ah 

1998 

Kobulnicky 

& Koo 

2000; Mehlert et ah 2002J Lemoinc-Busserollc et ah 

2003 

Erb et ah 2006; 


Maier et ah 2006). 


Here, we aim to discuss the metal enrichment histories of galactic systems in 
a hierarchical clustering scenario. The results of the simulations, the comparison 


with the data, and the implications of our results are described in Sect. |5.2[ We 
summarise our conclusions in Sect. 


5.2 Results 


We have constructed an ensemble of 112 simulated late-type galaxies in a chemo- 
dynamical framework (Section |3.2[ ), spanning a factor of ~ 50 in total mass, and 
sampling a range of assembly histories at a given total mass. 

This sample includes the same 89 semi-cosmological GCD+ simulations analysed 
in Renda et ah (2005b) and 23 more simulations, completed within the same frame¬ 


work, with the same parameters described in Section T2 and distributed among 
the same mass “bins”. These latter simulations sample more patterns of small-scale 
density fluctuations, which lead to different hierarchical assembly histories. A sub¬ 
set of 9 simulations (with M tot = lO n M 0 , collapse redshift z c = 1.5 and spin 
parameter A = 0.054) samples more density fluctuations patterns. 

As a benchmark for the semi-cosmological framework, we have also analysed 


2 among the disc galaxy simulations presented in Bailin et ah (2005), which are 


fully cosmological GCD+ runs and adopt a A-dominated CDM cosmology (Bailin et 


ah 2005 and references therein). The photometric properties of a stellar particle in a 
simulation are modeled as those of a Simple Stellar Population (SSP hereafter) with 
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the same age and metallicity. The evolution of the SSP properties as a function of 


age and metallicity is taken from Mouhcine & Langon (2003). Optical properties 


have been converted into the SDSS photometric system as in Fukugita et al. (1996) 


In the following, the stellar mass is measured within the 15 kpc central region, 
at each redshift. Metallicity and age of the integrated stellar populations, and the 
gas-phase metallicity, are all measured within the 10 kpc central region, at each 
redshift. 

We estimate the Oxygen abundance as 12 + log(0/H). Caution should be taken 
when comparing the estimated abundances for simulated galaxies with the observed, 


since the latter, based on integrated spectra, may suffer from systematics (e.g., Mous- 


takas & Kennicutt 2006). 


5.2.1 Gas—phase Oxygen abundance vs. stellar mass 


Fig. 5T shows a correlation between gas-phase Oxygen abundance and stellar mass 
in the simulations, at different redshifts. The consistency between semi-cosmological 


and cosmological simulations (see also Section 4.2) is reassuring as a robustness test 
for the semi-cosmological framework. 

Massive galaxies with comparable stellar masses span a range of gas-phase abun¬ 
dances, due to their differing assembly histories, thus suggesting that stellar mass 
should not be considered a good tracer of gas-phase metallicity. 

The results at redshift z = 0 are compared and in broad agreement with the ob¬ 
served relation between stellar mass and gas-phase Oxygen abundance for a sample 
of SDSS galaxies at z ~ 0.1 (Tremonti et al. 2004). The dispersion of the gas-phase 


metallicity at a given stellar mass at z = 0 is similar to what is locally observed in 
the SDSS. This is a remarkable result, given that simulations and observations adopt 
different pipelines to estimate the gas-phase metallicity: for example, the SDSS fi¬ 
bre spectra preferentially sample the inner ~ 25% of a galaxy (e.g., 


Tremonti et 


al. 2004 and references therein) whereas we are considering the Oxygen abundance 


within the inner 10 kpc in the simulations. 

The collapse-redshift z c = 1.5 runs at M to t = 10 11 M© are broadly as 

evolved as the most metal-poor z c = 2.0 simulations at the same M tot . If the 
semi-cosmological simulations were constrained to agree with the observations at 
redshift z ~ 0, the collapse-redshift range at M tot = 10 11 M© would then be 


1.5 < 


< 2.0, thus suggesting that the observed low stellar mass galaxies may 
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Figure 5.1: Gas-phase Oxygen abundance against integrated stellar mass in the 
simulations, at different redshifts. The semi-cosmological simulations with collapse 
redshift z c = 2 are shown as filled symbols (heptagons for M tot = 5 x 10 12 M 0 , 
hexagons for M tot = 10 12 M 0 , pentagons for M tot = 5 x 10 11 M 0 , boxes for 
M tot = 10 11 M 0 ). Larger empty boxes for collapse redshift z c = 1.5 runs. The 
cosmological simulations are shown as large empty triangles. The z = 0 gas-phase 
Oxygen abundance in the simulations is compared with z ~ 0.1 galaxies in the SDSS 


(Trcmonti et ah, 2004): solid line for the median; dotted line for 16 and 84 percentile; 


dashed line for 2.5 and 97.5 percentile. The results at low and intermediate redshifts 


are compared with the samples observed in Savaglio et al. (2005) and Liang et al 


(2006), shown as five and six vertices stars, respectively. 
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Figure 5.2: Gas-phase Oxygen abundance against B-band luminosity. Symbols are 
the same as in Fig. 5.1 The redshift evolution in the simulations is compared with 


the relationship between the observed luminosity and the gas-phase metallicity at 


intermediate redshifts in Lilly et al. (2003), Liang et al. (2004), and Mouhcine et al. 
(2006), shown as four, five and six vertices stars, respectively. 
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Figure 5.3: The V-band luminosity-weighted metallicity of the integrated stellar 
populations against the integrated stellar mass, at different cosmic epochs. Symbols 
are the same as in Fig. 5.1. The results at z = 0 are compared with the z ~ 0.1 re¬ 


lation between stellar mass and metallicity in the SDSS (Gallazzi et al. 2005): solid 
line for the median; dotted line for 16 and 84 percentile. 
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Figure 5.4: The V-band luminosity-weighted age of the integrated stellar popula¬ 
tions, against the stellar mass at redshift z = 0. Symbols are the same as in Fig. 13 
The results are compared with the z ~ 0.1 relation between stellar mass and age 
in the SDSS (Gallazzi et al. 2005): solid line for the median; dotted line for 16 and 
84 percentile. 
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Figure 5.5: The z = 0 Colour-Magnitude relation in the simulations, 
are the same as in Fig. |5.1 


local CMR for a sample of disc-dominated SDSS galaxies selected by [Piz agno et al. 
(2005), shown as empty boxes. 
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have collapsed out of density peaks with smaller values of 5; (see also Section 3.2). 


Fig. 5.1 shows that the redshift evolution of the mass-metallicity relationship 
depends on galaxy mass. Gas-phase abundances do not significantly evolve in simu¬ 
lations with stellar mass larger than ~ 10 10 M 0 , and the mass-metallicity relation¬ 
ship remains nearly as flat as it is at z ~ 1. On the contrary, there is a systematic 
metal enrichment for low stellar mass simulations, which increase their stellar con¬ 
tents by a large factor, migrating from stellar masses of 10 9 M 0 to ~ 10 10 M 0 . 
Such behaviour is consistent with a scenario where massive galaxies have burnt a 
significant fraction of their gas reservoirs and are dominated by stars early on, thus 
leaving their metal enrichment almost unchanged since z ~ 1. Low stellar mass 
galaxies would rather form stars in a more extended fashion, and a sizeable fraction 
of their baryons would still be gaseous at the present epoch. 

For intermediate redshifts, the mass-metallicity relation in the simulations is 


compared and in broad agreement with that reported at 0.4 < z < 1 in Savaglio et 


al. (2005) and Liang et al. (2006), which too support a flat mass-metallicity relation 
for massive galaxies since z ~ 1 down to z ~ 0. 

The gas-phase Oxygen abundances for low stellar masses in the simulations agree 
with the observed values, whereas the metallicity for more massive simulations tend 
to be systematically higher than the observed abundances. Part of the offset could 
be explained by systematics affecting the measurements, especially at the high- 
metallicity end (e.g.: Ellison & Kewley 2005 Kennicutt et al. 2003). 


5.2.2 Gas—phase Oxygen abundance vs. galaxy luminosity 


As a reflection of the relationship between stellar mass and gas-phase metallicity, the 


relationship between gas metallicity and galaxy luminosity is displayed in Fig. 5.2 


The luminosity-metallicity relationship is compared and in broad agreement with 


observations at intermediate redshift: Lilly et al. (2003) at 0.5 < z < 0.9; Liang et 


al. (2004) at 0 < z < 1.2; Mouhcine et al. (2006) at 0.2 < z < 0.8. As previously 


noted, this is a remarkable result, given that simulations and observations adopt 
different pipelines to estimate the gas-phase metallicity. 

The Oxygen abundances of luminous galaxies in the simulations, i.e. B-band 
magnitude brighter than -20, are larger than what is observed at 0.4 < z < 0.7. 


Systematics in the observations may play a role here too, as discussed above in 
connection with the mass-metallicity relationship for massive galaxies. 
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5.2.3 Stellar metallicity/age vs. stellar mass 


Age and metallicity of integrated stellar populations are powerful tracers of galaxy 
star formation and metal enrichment histories. Fig. 5.3| shows the relationship be¬ 
tween stellar mass and V-band luminosity-weighted metallicity, at different red- 
shifts. 

The general behaviour is similar to that of the gas-phase Oxygen abundance as 
a function of stellar mass (Fig. |5.f ), although stellar metallicities are less dispersed 
than gas-phase Oxygen abundances at a given stellar mass. 

The results at redshift z = 0 in the simulations are compared with the observed 
relation for a large sample of SDSS galaxies at 0.005 < z < 0.22 as reported in 


Gallazzi et al. (2005), whose sample includes both quiescent early-type and actively 


star-forming galaxies, which might explain the smaller stellar metallicity disper¬ 
sion at a given stellar mass in the simulations (which have late-type morphologies) 
compared with the dispersion observed in the more heterogeneous SDSS sample. 


Fig. |5.4| shows the z = 0 relationship between stellar mass and V-band luminosity 
weighted age of the integrated stellar populations in the simulations, compared with 
the observed relationship. We find no relation between age and stellar mass. Massive 
galaxies with stellar mass larger than ~ 10 10 M 0 have age ranging from ~ 6.5 Gyr 
to ~ 9 Gyr, in broad agreement with the age derived for SDSS galaxies with com¬ 
parable stellar masses. For less massive galaxies in the simulations, the age range 
is broader, going from ~ 3 to ~ 10 Gyr, thus older than in Gallazzi et al. (2005), 


whose sample is however much less populated and thus less robust at stellar masses 
M* < 10 10 M 0 (see Fig. 6 in Gallazzi et al. 2005). 


The age issue in the simulations as outlined above signals that feedback should 
more effectively prevent gas from collapsing and cooling in the early stages of the 
hierarchy in order to ensure that most of the stellar content in low stellar mass 
galaxies is not formed early on. This has been already invoked in previous studies 


(e.g.: Steinmetz & Navarro 1999). 


5.2.4 The Colour—Magnitude Relation 


Another signature of feedback-related problems in our sample, Fig. 5.5| shows the 
z = 0 relationship between the (g - r) colour and the i-band luminosity in the 
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simulations, compare^] with a sample of disc-dominated local galaxies in the SDSS, 
selected by Pizagno et al. (2005). 

The integrated stellar population colours of luminous galaxies in the simulations 
are in good agreement with the observed colours at comparable luminosities, which is 
expected due to the agreement, in massive simulations, of age and stellar metallicity 
with the observations. The colour dispersion in the simulations - which reflects 
the diverse formation histories which have been sampled - agrees with the observed 
dispersion. This is a remarkable result, given that the observed sample has explicitly 
been selected to include disc-dominated galaxies with a disc-to-total luminosity 
ratio fa > 0.9 (see Pizagno et al. 2005), whereas the simulated sample consists of 
late-type galaxy simulations which have not been morphologically selected. 

Although the colour dispersion for less massive simulations is comparable with 
the observed dispersion, fainter simulations show integrated (g - r) colours generally 
redder than the observed, i.e. the simulated CMR is shallower than the observed. 
The discrepancy between observed and simulated colours for faint galaxies is due 
to the age discrepancy for faint galaxies between observations and simulations. Al¬ 
though faint simulations in our sample agree with stellar and gaseous metallicities 
for locally observed galaxies with comparable luminosities and/or stellar masses, the 
former are older, resulting in redder colours. 


5.3 Conclusions 


We have presented the analysis of a sample of semi-cosmological late-type galaxy 
simulations with stellar mass larger than ~ 10 9 M 0 , to explore galaxy metal enrich¬ 
ment over the last half of the age of the Universe. A particular emphasis has been 
placed upon the relationship between stellar mass and metallicity of both gaseous 
and stellar components. 

The main results can be summarised as follows. 


- A correlation between gas-phase Oxygen abundance and stellar mass is present 
at all redshifts since z ~ 1 down to z = 0. 


L Since dust-correction has not been implemented in our code, we have compared the colour- 


magnitude relation (CMR hereafter) with the sample of Pizagno et al. (2005), whose internal dust 


reddening was corrected for as in Tully et al. (1998), rather than comparing the CMR in the 


simulations with the low redshift New York University Value-Added Galaxy Catalogue of 


et al. 


et al. 


(2005), which do not take into account extinction but the Galactic extinction as in 


(1998). 


Blanton 


Schlegel 
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The z = 0 relationship between stellar mass and gas-phase metallicity in the 
simulations is in broad agreement with the z ~ 0.1 correlation observed in 
SDSS galaxies. 

The z = 0 dispersion of the gas-phase Oxygen abundances for simulations 
with comparable masses and diverse merging histories is in broad agreement 
with the z ~ 0.1 dispersion observed in SDSS galaxies. This suggests that 
the latter could be due to the differing formation histories of galaxies with 
comparable masses. 

The mass-metallicity relationship for massive simulations is ~ flat since z ~ 1 
down to z = 0, thus suggesting that they experience substantial metal enrich¬ 
ment early on. 

No relationship is found in the simulations between age and stellar mass. The 
age range for massive simulations is similar to that observed for comparable 
mass z ~ 0 galaxies. However, the z = 0 stellar populations of less massive 
simulations tend to be dominated by stars older than the observed. A mecha¬ 
nism to delay star formation seems necessary to solve the discrepancy, as also 
suggested by the z = 0 Colour-Magnitude Relation in the simulations, which 
is shallower than the observed. 



Chapter 6 


Conclusions and Future Directions 


Which is “normal”?...Is the Milky Way an outlier? Is there any stellar halo-galaxy 
formation symbiosis? These driving questions in mind, we have tried to face the 
challenge they pose and undertaken the path which has led to this PhD Thesis. Here 
we summarise what (we think that) we see from our point of view. 


We have presented in Chapter [3] (which synthesises the background which is in 
Appendices |B| - [H| see also jRenda et al. 2005b) an analysis of several late-type 
galaxy stellar halos within a grid of semi-cosmological simulations, with particular 
emphasis placed upon the relationship between stellar halo metallicity and galactic 
luminosity. It helps to stress here that “Halo Semantics” is currently a contro¬ 
versial topic (e.g., Ibata et al. 2007 and references therein): we have labelled as 
“topographical halo” the ensemble of stellar particles in a simulation at a projected 
radius R > 15 kpc, at z = 0. Although future analysis of simulations at resolutions 
higher than the current would further improve the study we have undertaken, we 
have shown that - at any given galaxy luminosity or, conversely, galaxy dynamical 
mass - halo metallicities in the simulations span a range in excess of ~ 1 dex, a 
result which is strengthened by the robustness tests we have performed. 


We suggest that the underlying driver of the halo metallicity dispersion can be 
traced to the diversity of galactic mass assembly histories inherent within the hi¬ 
erarchical clustering paradigm. Galaxies with a more protracted assembly history 
possess more metal-rich and younger stellar halos, with an associated greater dis¬ 
persion in age, than galaxies which experience more of a monolithic collapse. 

For a given galaxy luminosity (or dynamical mass), those galaxies with more 
extended assembly histories also possess more massive stellar halos, which in turn 
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leads to a direct correlation between the stellar halo metallicity and its surface 
brightness (as anticipated by earlier semi-analytical models - e.g.: Chapter [2 Renda 


et ah 2005). By extension, such a correlation may prove to be a useful diagnostic 


tool for disentangling the formation history of late-type galaxies. 


Recently, Mouhcine et al. (2005) have presented an observed correlation between 
stellar halo metallicity and galactic luminosity. The observed halo metallicity disper¬ 
sion at a given galactic luminosity is smaller than what we find in our simulations. 
However the latter can account for the outliers in the observed trend. Since our 
motivation has been to study which is the effect of the pattern of the initial density 
fluctuations alone on the stellar halo features at redshift z = 0 in late-type galaxy 
simulations, it helps to note that galaxy formation, as it is observed, is an ongoing 
process which is the result of the interplay among different parameters, of which the 
pattern of the initial density fluctuations (thus the merging history) is one. We have 
shown that the merging history alone may be held responsible of the dispersion in 
halo metallicity at comparable galactic luminosities. 

We have then analysed the stellar mass assembly, and the relationships between 
stellar, gaseous and total mass since z ~ 1 down to z = 0 in the simulations, to 
explore whether hierarchical formation scenarios may account for the constraints on 
galaxy mass assembly over the second half of the age of the Universe (Chapter [4]). 

The redshift evolution of the stellar mass in simulations with comparable lumi¬ 
nosities at z = 0 shows a large dispersion, which lies in their differing merging 

histories. This contrasts with the redshift evolution of the total mass within the 
central 100 kpc, which is already assembled by z ~ 1.5. The stcllar-to-total mass 


ratio in the simulations is in broad agreement with the observed (e.g., Conselice et 


al. 2005). We find that massive simulations are generally more evolved than their 


lower mass counterparts. Galaxies with comparable luminosities/stellar masses at 
z = 0 in the simulations can have significantly different assembly histories of their 
stellar content, and diverse gaseous contents. This suggests that stellar mass should 
not be considered neither a robust tracer of the galaxy merging history nor a good 
tracer of the total baryonic galaxy mass. 

We have finally analysed the metal enrichment of the stellar and gaseous com¬ 
ponents in the simulations (Chapter [5]). A correlation between gas-phase Oxygen 
abundance and stellar mass is present since z ~ 1 down to z = 0. The z = 0 re¬ 
lationship between stellar mass and gas-phase metallicity in the simulations is in 
broad agreement with the z ~ 0.1 correlation observed in SDSS galaxies (Tremonti 















et al. 2004). The z = 0 dispersion of the gas-phase Oxygen abundances for simu¬ 


lations with comparable masses and diverse merging histories is in broad agreement 
with the dispersion observed at z ~ 0.1 in SDSS galaxies. This suggests that the 
latter could be due to the differing formation histories of galaxies with comparable 
masses. The mass-metallicity relationship for massive simulations is ~ flat since 
z ~ 1 down to z = 0, thus suggesting that they experience substantial metal enrich¬ 
ment early on. No relationship is found in the simulations between age and stellar 
mass. The age range for massive simulations is similar to that observed for compa¬ 
rable mass z ~ 0 galaxies. However, the z = 0 stellar populations of less massive 
simulations tend to be dominated by stars older than the observed. A mechanism to 
delay star formation seems necessary to solve the discrepancy, as also suggested by 
the z = 0 Colour-Magnitude Relation in the simulations, which is shallower than 
the observed. 

Where do we go from here? However promising they are, chemo-dynamical 
simulations still suffer from feedback related issues (e.g.: Steinmetz & Navarro; j j l999[ 
Thacker & Couchmanl 120001 lAbadi et al. 1 120031 IRobertson et al. 1 120061 IStinson et 


al. 2006 Kaufmann et al. 2007) which are reflected here in the mass-metallicity 


relationship and in the Colour-Magnitude Relation (Chapter [5]). It is reassuring 
though - as a robustness test for the semi-cosmological framework - that there is 
consistency between the results of semi-cosmological and cosmological simulations, 
as shown in Chapters [| and [5j 

We have focused here in this PhD Thesis on the formation of stellar halos in 
late-type galaxies. We have shown that the diverse galaxy merging histories may be 
held responsible of the observed dispersion in halo metallicity at comparable galactic 
luminosities (Chapter [3]). Higher resolution simulations, a deeper understanding of 
the role of feedback processes in galaxy formation, and further observations are all 
together needed to take the next step on this research path. 
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Appendix A 

On the origin of fluorine 
in the Milky Way 


The main astrophysical factories of fluorine ( 19 F) are thought to be Type II super¬ 
novae, Wolf-Rayet stars, and the Asymptotic Giant Branch (AGB) of intermediate 
mass stars. We present a model for the chemical evolution of fluorine in the Milky 
Way in a semi-analytic multi-zone chemical evolution framework. For the first time, 
we quantitatively demonstrate the impact of fluorine nucleosynthesis in Wolf-Rayet 
and AGB stars. The inclusion of these latter two fluorine production sites provides 
a possible solution to the long-standing discrepancy between model predictions and 
the fluorine abundances observed in Milky Way giants. Finally, fluorine is discussed 
as a possible probe of the role of supernovae and intermediate mass stars in the 
chemical evolution history of the globular cluster oo Centauri. 


A.l Introduction 


The three primary astrophysical factories for fluorine ( 19 F) production have long 
been thought to be Type II Supernovae (SNe II), Wolf-Rayet (WR) stars, and 
Asymptotic Giant Branch (AGB) stars (e.g., respectively: Woosley & Weaver||1995{ 


Meynet & Arnould 2000; Forestini et al. 1992 Mowlavi, Jorissen & Arnould 1998). 
Previous attempts to model the Galactic production and evolution of 19 F have been 
restricted to explore the role of SNe II alone (e.g.: Timmes, Woosley & Weaver 1995; 
Alibes, Labay & Canal 2001). 

The above problem has now been ameliorated by the release of the first detailed 
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yield predictions for fluorine production from WR and AGB stars. We are now in a 
position to incorporate these yields into a Galactic chemical evolution framework, in 
order to assess the respective contributions of the three putative fluorine production 
sites. To do so, we will make use of GEtool, a semi-analytical multi-zone Galactic 
chemical evolution package which has been calibrated with extant observational data 


for the Milky Way ( Fenner fe Gibson] 2003 Gibson et al. 2003). 

Specifically, in what follows, we compare the model fluorine distribution in the 
Milky Way with the abundances observed by Jorissen, Smith & Lambert (1992) in 
near-solar metallicity giants. Further, our model predictions are contrasted with 
new fluorine determinations for giants in the Large Magellanic Cloud (LMC) and 


uj Centauri 

Cunha et al. 

2003 

). In addition, new results for more u ; Centauri 

giants from 

Smith et al. ( 

2004 

) are included. The latter two systems are likely to 


have had very different star formation and chemical evolution histories from those 
of the Milky Way, but despite these obvious differences, a comparison against these 
new data can be valuable. In Section A.2 we provide a cursory overview of the 
three traditional 19 F nucleosynthesis sites; the chemical evolution code in which the 
nucleosynthesis products from these factories have been implemented is described 
in Section [A. 3 Our results are then presented and summarised in Sections A.4 and 
A.5 respectively. 


A.2 Nucleosynthesis of 19 F 


A.2.1 Type II Supernovae 

The massive star progenitors to SNe II produce fluorine primarily as the result 


of spallation of 20 Ne by /i and r neutrinos near the collapsed core (Woosley & 


Haxton 1988 Woosley et al. 1990). A fraction of the 19 F thus created is destroyed 


by the subsequent shock but most is returned to the ambient Interstellar Medium 
(ISM). The fluorine yields by neutrino spallation are very sensitive to the assumed 


spectra of jd and r neutrinos (Woosley et ah. 2002), which could be nonthermal 


and deficient on their high-energy tails, lowering the equivalent temperature of the 


neutrinos in the supernova model (Myra & Burrows, 1990). An additional source 
of 19 F derives from pre-explosive CNO burning in helium shell. However, fluorine 
production by neutrino spallation is largely dominant, as evident by comparing 


the models in Woosley & Weaver (1995), and recent models which do not include 
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neutrino nucleosynthesis of fluorine (Limongi & Chieffi 2003). Most recently, Heger 


et al. (2004) suggest that the relevant neutrino cross sections need to be revised 
downwards; if confirmed, the associated SNe 11 19 F yield would decrease by ~ 50%. 
In light of the preliminary nature of the Heger et al. claim, we retain the conservative 


choice offered by the Woosley & Weaver (1995) compilation. 


A.2.2 Asymptotic Giant Branch stars 

The nucleosynthesis pathways for fluorine production within AGB stars involve both 


helium burning and combined hydrogen-helium burning phases (e.g.: Forestini et al 


1992 Jorissen et al.|1992 Mowlavi et al.|1998) and are companions for the nucleosyn¬ 


thesis by slow neutron accretion (s-process) (Mowlavi et al. 1998). Provided a suit¬ 
able source of protons is available, fluorine can be synthesised via 14 N(a, 7 ) 18 F (/3 + ) 
18 0(p,a) 15 N(a, 7 ) 19 F. Primary sources of uncertainty in predicting fluorine nucle¬ 
osynthesis in AGB stars relate to the adopted reaction rates, especially 14 C(ct, 7) 18 0 
and 19 F(cp p) 22 Ne, and the treatment of the nucleosynthesis occurring during the 
convective thermal pulses. Nucleosynthesis during the interpulse periods can also 
be important if protons from the envelope are partially mixed in the top layers 


of the He intershell (partial mixing zone), as Lugaro et al. (2004) have recently 


demonstrated. Nucleosynthesis in this zone may result in a significant increase in 
the predicted 19 F yields. The magnitude of these systematic uncertainties for stellar 
models with mass ~ 3 M 0 and metallicities Z = 0.004 - 0.02 are ~ 50 %, while 
for stellar models with mass M = 5 M 0 and metallicity Z = 0.02 the uncertainty 
is a factor of ~ 5, due to the uncertain 19 F(a, p) 22 Ne reaction rate. Characterising 
the mass- and metallicity- dependence of the partial mixing zone- 19 F relationship 
needs to be completed before we can assess its behaviour self-consistently within 
our chemical evolution model of the Milky Way. For the present study, we have 


adopted the yields presented in Renda et al. (2004), based upon the Karakas & 


Lattanzio (2003, and references therein) models, which themselves do not include 
19 F nucleosynthesis via partial mixing. This choice is a conservative one, and thus 
should be considered as a lower limit to the production of 19 F from AGB stars. 

For stars more massive than « 4 M 0 , the convective envelope is so deep that it 
penetrates into the top of the hydrogen-burning shell so that nucleosynthesis actu¬ 
ally occurs in the envelope of the star. Such “hot-bottom-burning” acts to destroy 
19 F, and should be treated self-consistently within the AGB models considered. 
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A.2.3 Wolf—Rayet stars 


Fluorine production in WR stars is tied to its nucleosynthesis during the helium¬ 
burning phase. At the end of this phase though, significant fluorine destruction 
occurs via 19 F(cq p) 22 Ne. Any earlier synthesised 19 F must be removed from the 
stellar interior in order to avoid destruction. For massive stars to be significant con¬ 
tributors to net fluorine production, they must experience mass loss on a timescale 
that allows the removal of 19 F before its destruction. This requirement is met by 
WR stars. 


Recently, Meynet & Arnould (2000) studied the role that such stars can play in 


the chemical evolution of fluorine by adopting updated reaction rates coupled with 
extreme mass-loss rates in not-rotating stellar models. They pointed out that WR 
mass-loss is strongly metallicity-dependent, and that the number of WR stars at 
low metallicities is very small. Their WR yields reflect such metallicity-dependence, 
with minimal fluorine returned to the ISM at low metallicities, but significant 19 F 
returned at solar and super-solar metallicities. The WR yields are sensitive to 
the adopted reaction and mass-loss rates, while rotating models could favour an 
early entrance into the WR phase for a given mass, decrease the minimum initial 
mass for a star to go through a WR phase at a given metallicity, and open more 
nucleosynthetic channels because of the mixing induced by rotation. Therefore, 


after Meynet & Arnould (2000), we consider the aforementioned WR yields as lower 
limits. 


A.3 The model 


In this study we employ GEtool, our semi-analytical multi-zone chemical evolution 


package to model a sample Milky Way-like disk galaxy (Fenner & Gibson 2003 


Gibson et al. 2003). A dual-infall framework is constructed in which the first infall 


episode corresponds to the formation of the halo, and the second to the inside-out 
formation of the disk. 

A Kroupa, Tout & Gilmore (1993) initial mass function (IMF) has been assumed, 
with lower and upper mass limits of 0.08 M 0 and 120 M e , respectively. Stellar yields 
are one of the most important features in galactic chemical evolution models, yet 
questions remain concerning the precise composition of stellar ejecta, due to the 
uncertain role played by processes including mass loss, rotation, fall-back, and the 
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location of the mass cut, which separates the remnant from the ejected material in 


SNe. The SNe II yields are from Woosley & Weaver (1995); for stars more massive 


than 60 M 0 , the yields are assumed to be mass-independent. Such assumption is 
made to avoid extreme extrapolation from the most massive star in the Woosley & 
Weaver models (40 M 0 ) to the upper end of the IMF (120 M 0 ), and has negligible 
effect on the results, given the shape of the adopted IMF. 


We have halved the iron yields shown in Woosley & Weaver (1995), as suggested 


by Timmes et al. (1995). The Type la (SNe la) yields of Iwamoto et al. (1999) 


were also employed. We adopted the metallicity-dependent yields of Renzini & Voli 
(1981) for single stars in the mass range 1-8 M 0 . For the purposes of this work, 


which focuses on fluorine, the choice of the Renzini & Voli yields set does not affect 


the results. Metalllicity-dependent stellar lifetimes have been taken from Schaller 


et ah (1992) 


We have constructed three Milky Way (MW) model variants that differ only in 
their respective treatments of 19 F production: 1) MWa assumes that SNe II are the 
only source of 19 F; 2) MWb includes yields from both SNe II and WR stars; 3) MWc 
includes all three sources of fluorine - SNe II, WR, and AGB stars. 

We end by noting that within our adopted dual-infall framework for the chemical 
evolution of the Milky Way, our model is constrained by an array of observational 
boundary conditions, including the present-day star and gas distributions (both in 
density and metallicity), abundance ratio patterns, age-metallicity relation, and G- 


dwarf distribution (Gibson et ah, 2003). While the modification of any individual 


ingredient within the model framework will have an impact, to some degree, upon 
the predicted chemical evolution, this can only eventuate at the expense of one 
or more of the aforementioned boundary conditions that we require our model to 
adhere to. Within our framework, yield uncertainties will dominate the systematic 
uncertainties for the predicted evolution of 19 F. 


A.3.1 Fluorine yields 

We now summarise the 19 F yields employed in our three “Milky Way” models. 


1) SNe II 19 F yields are taken from Woosley & Weaver (1995) and assumed to 


be mass-independent for stellar masses in excess of 60 M 0 . 


2) WR 19 F yields are taken from Meynet & Arnould (2000) for stellar masses in 
the range 25 - 120 M 0 : each star within this range is assumed to evolve through the 
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WR stage. Such simplifying assumption could overestimate the WR contribution 
to fluorine, even though the adopted WR yields are themselves lower limits (Sec¬ 
tion A.2.3). The WR fluorine contribution has been added to the corresponding 
SNe II contribution (which comes from a different stage of the stellar evolution). 

3) AGB 19 F and oxygen yields in the mass range 1 - 6.5 M 0 have been de¬ 
rived from stellar models constructed with the Mount Stromlo Stellar Structure 


Code (Frost & Lattanzio 1996; Karakas et al. 2002), and are presented in Renda et 


al. (2004). The post-processing nucleosynthesis models with 74 species and time- 


dependent diffusive convective mixing are described in detail in Frost et ah (1998) 


and Karakas & Lattanzio (2003). 

To ensure internal consistency, we have also employed the AGB oxygen yields in 
lieu of those of Renzini & Voli (1981), within this mass range. 

The above fluorine yields are shown in Figure A.l In Figure A.2 the yields are 
expressed as [F/Ojj]and ([F/0 ])/mf, the latter corresponding to the mean [F/O] 
yields for SNe II and AGB stars, weighted by the IMF over the SNe II and AGB 
mass range, respectively. We have not shown a comparable entry for the WR stars 
as a self-consistent treatment of the oxygen production was not included in Meynet 
& Arnould (2000). Here, oxygen has been used as the normalisation to make easier 
the comparison with the observations, especially in c o Centauri, though oxygen can 
be synthesised in various stellar sites, and its yields can be affected by different 
reaction rates and modeling of helium cores, semi-convection, convective boundary 
layers, and mass-loss (e.g.: Woosley et ah 2002 Dray et ah 2003). 


A.4 Results 


In Figure A.3 the evolution of [F/O], A(O), the gas infall rate ctinfaiu the star 
formation rate (SFR), the SNe II rate and the gas-phase global metallicity Z of 
the three models at the solar neighbourhood are summarised. The empirical SFR 
history derived by Bertclli & Nasi (2001) is shown as a thick solid line in panel “d” 
of Figure |iO] while the shaded region corresponds to the range of values suggested 
by Rana (1991). A conservative range of estimated SNe II rates is also shown 


1 Hereafter, [X/Y] = log 10 (X/Y) — log lo (X/Y) 0 and A(X) = 12 + log 10 (nx/iiH). An accurate 
determination of photospheric solar abundances requires detailed modeling of the solar granulation 
and accounting for departures from local thermodynamical equilibrium (e.g., Allende Prieto, Lam¬ 
bert & Asplund 2001). We adopt the solar fluorine abundance suggested by Cunha et al. (2003), 
and the solar iron and oxygen abundances from Holweger (2001). 
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Figure A. 2: [F/0] and ([F/0])imf for SNe II and AGB yields (upper and lower 
panels, respectively). Here, A( 19 F) 0 = 4.55 (see discussion in Cunha et al. 2003) 
and A(O) 0 = 8.736 (e.g., Holweger 2001). The shaded regions in panels “a” and 
“b” show the observed [F/O] in u Cen giants (Cunha et al. 2003). The ([F/0])imf 
are weighted by the IMF over the SNe II (11 - 40 M 0 ) and AGB (1 - 6.5 M e ) 
mass range, respectively. In panel “d”, both ([F/0]agb)imf and ([F/0]sNe ii)imf 
are shown (closed boxes and open triangles, respectively). 
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Figure A.3: Predicted evolution in the solar neighbourhood of: “a”, the gas infall 
rate a in f a if, “b”, A(O); “c”, [F/O]; “d”, the star formation rate (SFR); “e”, the 
SNe II rate; “f”, the metallicity Z (MWa, solid line; MWb, dotted; MWc, short- 
dashed). The SFR history at the solar neighbourhood obtained by Bertelli & Nasi 
(2001) is also shown as a thick solid line in panel “d”, while the shaded region shows 
the range of values suggested by Rana (1991). A range of values corresponding to 
the estimated SNe II rate is shown in panel “e” (Cappellaro et al. 1999). 
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in panel “e” (Cappcllaro, Evans & Turatto 1999) Q Figure A.4 then shows the the 
evolution of [F/0] versus A(0) (panel “a”), and the evolution of [F/0] versus [O/Fe] 


(panel “b”), compared against the IMF-weighted SNe II yields (recall Figure A.2). 

The MWa model provides a satisfactory reproduction of the estimated star for¬ 
mation history and SNe II rate in the solar neighbourhood (Figure A.3 panels “d” 
and “e”). This model, whose only fluorine source is SNe II, underproduces fluorine 
with respect to the abundances measured in K and M Milky Way giants observed by 


Jorissen et ah (1992) and reanalysed by Cunha et al. (2003) (Fig. A.4 panel “a”). 


Fig. A.4 panel “a” does not show the s-process enriched AGB stars of spectral 
types MS, S, or C in Jorissen et al. (1992), where freshly synthesised fluorine could 
be mixed to the stellar surface. Such inclusion of self-polluted 19 F-rich stars could 
obscure any metallicity trend. The results of the MWa and MWb models show that 
the additional contribution from WR stars increases [F/O] by up to factor of 2 by 
the present-day, but it is negligible in excess of ~ 9 Gyr ago (Figure [A~3| panel “c”). 

The addition of both WR and AGB sources within the MWc model leads to a 
present-day [F/O] that is ~ 0.4 dex greater than in the MWa case. Further, and 
perhaps more important, AGB stars are now shown to deliver significant amounts 
of fluorine to the ISM during the early epochs of the Milky Way’s evolution. Such a 
result is entirely consistent (and expected) given the metallicity-dependence of the 
AGB yields; said yields possess [F/O] ratios which are greater at lower metallicities 
(recall Figure [AY] ). We can conclude that it is only the addition of both the WR and 
the AGB contributions which allow for a significant improvement in the comparison 
between galactic models incorporating fluorine evolution and the observational data. 


A.5 Discussion 


We have studied the Galactic chemical evolution of fluorine, for the first time us¬ 
ing new grids of stellar models which provide self-consistent predictions of fluorine 
nucleosynthesis for stars in both the WR and AGB phases of stellar evolution. We 
have shown that the WR contribution is significant at solar and super-solar metallic¬ 
ities because of the adopted metallicity-dependent mass-loss prescription employed 


2 The range of values shown in panel “e” of Figure A.3 is derived from the sample of SOa - Sb 


galaxies in Cappellaro et al. (1999) - 0.42 ±0.19 SNu, where 1 SNu= 1 SN(100 yr) -1 (10 10 Lg) -1 , 
assuming Lf IW = 2 x 10 10 Lg and a galactic radial extent of 15 kpc. Given these assumptions, 
the estimated SNe II rate at the solar neighbourhood is necessarily uncertain. 
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Figure A.4: In panel “a”, [F/O] as a function of A(O) for the MW models (MWa, 
solid line; MWb, dotted; MWc, short-dashed). Also shown are the values observed 
in Milky Way, LMC, and c o Cen giants (crosses, boxes and hexagons, respectively). 
In panel “b”: [F/O] as a function of [O/Fe], compared with the IMF-weighted 
([F/0]gNe ii )imf an d ([0/Fe]sNe ii)imf yields for SNe II (open triangles). Within the 
open triangles: “0” corresponds to Z=0; “1”, to Z=1.9xl0“ 6 ; “2”, to Z=1.9xl0" 4 ; 
“3”, to Z=1.9xl0~ 3 ; “4”, to Z=1.9xl0 -2 . The upper panels represent enlargements 
of the framed regions delineated in the corresponding bottom panels. 
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in the stellar models. In contrast, the contribution of AGB stars to fluorine pro¬ 
duction peaks during the early epochs of the Galaxy’s evolution (again due to the 
metallicity-dependent behaviour of the AGB models). In combination, the addi¬ 
tion of the WR and AGB yields leads to a significant improvement in the galactic 
chemical evolution models when compared against observations. 

The comparison between our MW models and the fluorine abundances in LMC 
and oo Gen giants (Cunha et al. 2003) is not straightforward, as the latter two have 
star formation (and therefore chemical evolution) histories different from that of 
the MW. However, it is interesting to speculate on the possible origin of fluorine in 
oo Gen, given the unique nature of this “globular cluster” (e.g., Smith 2003). Specif¬ 
ically, co Gen is the most massive Galactic cluster, and unlike most globulars, pos¬ 
sesses a significant spread in metallicity (~ 1.5 dex) amongst its stellar population. 
It has been suggested that oo Cen is actually the remnant core of a tidally-disrupted 
dwarf galaxy (Bekkilh Freeman, 2003). Such a scenario could naturally drive radial 
gas inflows to the dwarf nucleus, potentially triggering starbursts. 


Interestingly, SNe II ejecta are characterised by low ([F/0 ])/mf (Figure A.2 
panel ”c”) and high ([0/Fe])/MF (Figure A.4| panel ”b”), whereas AGB ejecta have 
higher ([F/0 ])jmf (Figure A.2 panel “d”). The observed oo Cen giants have pri¬ 
marily low [F/O] (Figures: [W2| panels “a” - “b”; |A.4[ panels “a” - “b ”) and high 
[O/Fe] (Figure [A .4 , panel “b”) values, consistent with a picture in which their in¬ 
teriors have been polluted by the ejecta of an earlier generation of SNe II, but not 
from a comparable generation of AGB. Given that the observed oxygen abundance 
in such oo Cen giants is similar to that seen in comparable LMC and MW giants, as 
in Figure |A~4j panel “a”, this would suggest that the chemical enrichment of oo Cen 
proceeded on a short timescale (to avoid pollution from the lower mass progenitors 
to the AGB stars) and in an inhomogeneous manner (given the significant scatter 
in observed fluorine abundances), as previously discussed by Cunha et al. (2003). 
Should the (downward) revised neutrino cross sections alluded to in Section A.2.1 be 
confirmed (Heger et al. 2004), the concurrent factor of ~ 2 reduction in SNe II 19 F 
production would improve the agreement of the model with the observed [F/O] ratio 
in oo Cen giants. This would consequently strengthen our conclusions which already 
support a picture whereby these giants have been polluted by earlier generations of 
SNe II ejecta. 
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Figure B.l: MDFs for the halo fields observed in Mouhcine et ah (2005), reanalysed 
by the same pipeline as in Renda et ah (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF refers to are also shown. 
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Figure B.2: (continued) MDFs for the halo fields observed in Mouhcine et ah (2005), 
reanalysed by the same pipeline as in Renda et ah (2005b). The 68% Confi¬ 
dence Level range and the number of stellar particles each MDF refers to are also 
shown. 
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Figure B.3: Halo stellar particle MDFs for the M tot = 10 11 M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure B.4: Halo stellar particle MDFs for the M tot = 5xlO n M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure B.5: (continued) Halo stellar particle MDFs for the M to t = 5xlO n M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF refers to are also shown. 
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Figure B.6: (continued) Halo stellar particle MDFs for the M tot = 5xlO n M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF refers to are also shown. 
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Figure B.7: (continued) Halo stellar particle MDFs for the M to t = 5xlO n M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF refers to are also shown. 
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Figure B.8: Halo stellar particle MD1 
simulations in Renda et al. (2005b). 
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Figure B.9: (continued) Halo stellar particle MDFs for the M to t = 10 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF refers to are also shown. 
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Figure B.10: (continued) Halo stellar particle MDFs for the M tot = 10 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF refers to are also shown. 
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Figure B.ll: (continued) Halo stellar particle MDFs for the M tot = 10 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF refers to are also shown. 












Figure B.12: (continued) Halo stellar particle MDFs for the M toti = 10 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 


and the number of stellar particles each MDF refers to are also shown. 
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Figure B.13: Halo stellar particle MDFs for the M tot = 5xl0 12 M e semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF refers to are also shown. 
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Figure B.14: (continued) Halo stellar particle MDFs for the M tot = 5xl0 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF refers to are also shown. 







118 


APPENDIX B. STELLAR HALO METALLICITY DISTRIBUTIONS 



-4 -3 -2 -1 0 -4 -3 -2 -1 0 

[Fe/H] [Fe/H] 


Figure B.15: (continued) Halo stellar particle MDFs for the M to t = 5xl0 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF refers to are also shown. 
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Figure B.16: (continued) Halo stellar particle MDFs for the M tot = 5xl0 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF refers to are also shown. 







120 


APPENDIX B. STELLAR HALO METALLICITY DISTRIBUTIONS 



[Fe/H] [Fe/H] 


Figure B.17: (continued) Halo stellar particle MDFs for the M to t = 5xl0 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF refers to are also shown. 
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Figure C.l: Halo stellar particle [O/Fe] distribution for the M tot = 10 11 M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each distribution refers to are also shown. 
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Figure C.2: Halo stellar particle [O/Fe] distribution for the M tot = 5xlO n M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each distribution refers to are also shown. 
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Figure C.3: (continued) Halo stellar particle [O/Fe] distribution for the 
M tot = 5xlO n M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Confidence Level range and the number of stellar particles each distribution 
refers to are also shown. 
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Figure C.4: (continued) Halo stellar particle [O/Fe] distribution for the 
M to t = 5xlO n M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Confidence Level range and the number of stellar particles each distribution 
refers to are also shown. 
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Figure C.5: (continued) Halo stellar particle [O/Fe] distribution for the 
M tot = 5xlO n M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Confidence Level range and the number of stellar particles each distribution 
refers to are also shown. 





Figure C.6: Halo stellar particle [O/Fe] distribution for the M tot = lxlO 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each distribution refers to are also shown. 
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Figure C.7: (continued) Halo stellar particle [O/Fe] distribution for the 
M tot = 10 12 M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Confidence Level range and the number of stellar particles each distribution 
refers to are also shown. 
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Figure C.8: (continued) Halo stellar particle [O/Fe] distribution for the 
M to t = 10 12 M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Conhdence Level range and the number of stellar particles each distribution 
refers to are also shown. 
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Figure C.9: (continued) Halo stellar particle [O/Fe] distribution for the 
M tot = 10 12 M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Conhdence Level range and the number of stellar particles each distribution 
refers to are also shown. 
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Figure C.10: (continued) Halo stellar particle [O/Fe] distribution for the 
M to t = 10 12 M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Confidence Level range and the number of stellar particles each distribution 
refers to are also shown. 
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Figure C.ll: Halo stellar particle [O/Fe] distribution for the M tot = 5xl0 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each distribution refers to are also shown. 
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Figure C.12: (continued) Halo stellar particle [O/Fe] distribution for the 
M to t = 5xl0 12 M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Confidence Level range and the number of stellar particles each distribution 
refers to are also shown. 
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Figure C.13: (continued) Halo stellar particle [O/Fe] distribution for the 
M tot = 5xl0 12 M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Confidence Level range and the number of stellar particles each distribution 
refers to are also shown. 





Figure C.14: (continued) Halo stellar particle [O/Fe] distribution for the 
M to t = 5xl0 12 M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Confidence Level range and the number of stellar particles each distribution 
refers to are also shown. 
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Figure C.15: (continued) Halo stellar particle [O/Fe] distribution for the 
M tot = 5xl0 12 M 0 semi-cosmological simulations in Renda et al. (2005b). The 
68% Confidence Level range and the number of stellar particles each distribution 
refers to are also shown. 
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Figure D.l: Halo (red) and galaxy (black) age distributions for the M tot = 10 11 M Q 
semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.2: Halo (red) and galaxy (black) age distributions for the 

Mtot = 5xlO n M 0 semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.3: (continued) Halo (red) and galaxy (black) age distributions for the 
M to t = 5xlO n M 0 semi-cosmological simulations in Renda et al. (2005b). 


























145 



10 5 10 5 


age (Gyr) age (Gyr) 

Figure D.4: (continued) Halo (red) and galaxy (black) age distributions for the 
Mtot = 5xlO n M 0 semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.5: (continued) Halo (red) and galaxy (black) age distributions for the 
M to t = 5xlO n M 0 semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.6: Halo (red) and galaxy (black) age distributions for the M tot = 10 12 
semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.7: (continued) Halo (red) and galaxy (black) age distributions for the 
M tot = 10 12 M 0 semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.8: (continued) Halo (red) and galaxy (black) age distributions for the 
Mtot = 10 12 M 0 semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.9: (continued) Halo (red) and galaxy (black) age distributions for the 
M tot = 10 12 M 0 semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.10: (continued) Halo (red) and galaxy (black) age distributions for the 
Mtot = 10 12 M 0 semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.ll: Halo (red) and galaxy (black) age distributions for the 

M to t = 5xl0 12 Mq semi-cosmological simulations in Renda et al. (2005b). 

































153 



10 5 10 5 

age (Gyr) age (Gyr) 


Figure D.12: (continued) Halo (red) and galaxy (black) age distributions for the 
Mtot = 5xl0 12 M 0 semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.13: (continued) Halo (red) and galaxy (black) age distributions for the 
M to t = 5xl0 12 Mq semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.14: (continued) Halo (red) and galaxy (black) age distributions for the 
Mtot = 5xl0 12 M 0 semi-cosmological simulations in Renda et al. (2005b). 
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Figure D.15: (continued) Halo (red) and galaxy (black) age distributions for the 
M to t = 5xl0 12 Mq semi-cosmological simulations in Renda et al. (2005b). 
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Figure E.l: Galaxy rotation curves for the M tot = 5xlO n M 0 semi-cosmological 
simulations in Renda et al. (2005b). Black dots for the gaseous particles. Red dots 
for the stellar particles; red boxes for the ensemble of stellar particles (enclosed 
within the magenta frame) which are counter-rotating with vg < 0 within each 
topographical halo at a projected radius R > 15 kpc. The size of such ensemble is 
also shown in each panel. 
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Figure E.2: (continued) Galaxy rotation curves for the M tot = 5xlO n M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 
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Figure E.3: (continued) Galaxy rotation curves for the M tot = 5xlO n M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with vg < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 





























163 



R (kpc) R (kpc) 


Figure E.4: (continued) Galaxy rotation curves for the M tot = 5xlO n M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 
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Figure E.5: Galaxy rotation curves for the M tot = 10 12 M 0 semi-cosmological simu¬ 
lations in Renda et al. (2005b). Black dots for the gaseous particles. Red dots for the 
stellar particles; red boxes for the ensemble of stellar particles (enclosed within the 
magenta frame) which are counter-rotating with < 0 within each topographical 
halo at a projected radius R > 15 kpc. The size of such ensemble is also shown in 
each panel. 
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Figure E.6: (continued) Galaxy rotation curves for the M tot = 10 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 




























166 


APPENDIX E. GALAXY ROTATION CURVES 



0 20 40 60 80 

R (kpc) 


Figure E.7: (continued) Galaxy rotation curves for the M tot = 10 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with vg < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 
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Figure E.8: (continued) Galaxy rotation curves for the M tot = 10 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 
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Figure E.9: (continued) Galaxy rotation curves for the M tot = 10 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with vg < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 
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Figure E.10: Galaxy rotation curves for the M tot = 5xl0 12 M 0 semi-cosmological 
simulations in Renda et al. (2005b). Black dots for the gaseous particles. Red dots 
for the stellar particles; red boxes for the ensemble of stellar particles (enclosed 
within the magenta frame) which are counter-rotating with vg < 0 within each 
topographical halo at a projected radius R > 15 kpc. The size of such ensemble is 
also shown in each panel. 
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Figure E.ll: (continued) Galaxy rotation curves for the M tot = 5xl0 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with vg < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 
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Figure E.12: (continued) Galaxy rotation curves for the M tot = 5xl0 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 
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Figure E.13: (continued) Galaxy rotation curves for the M tot = 5xl0 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with vg < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 
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Figure E.14: (continued) Galaxy rotation curves for the M tot = 5xl0 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). Black dots for the gaseous parti¬ 
cles. Red dots for the stellar particles; red boxes for the ensemble of stellar particles 
(enclosed within the magenta frame) which are counter-rotating with < 0 

within each topographical halo at a projected radius R > 15 kpc. The size of such 
ensemble is also shown in each panel. 
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APPENDIX F. (ONCE MORE) STELLAR HALO MDF 
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Figure F.l: Kinematically selected halo MDFs for the M to t = 5xlO n M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF relates to are also shown. Each 
MDF refers to the stellar particles in the simulation which are counter-rotating 
with v 0 < 0 at a projected radius R > 15 kpc. 
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Figure F.2: (continued) Kinematically selected halo MDFs for the M tot = 5xlO n M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF relates to are also shown. Each 
MDF refers to the stellar particles in the simulation which are counter-rotating with 
v 0 < 0 at a projected radius R > 15 kpc. 
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Figure F.3: Kinematically selected halo MDFs for the M to t = 10 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF relates to are also shown. Each 
MDF refers to the stellar particles in the simulation which are counter-rotating 
with v 0 < 0 at a projected radius R > 15 kpc. 
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Figure F.4: (continued) Kinematically selected halo MDFs for the M tot = 10 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF relates to are also shown. Each 
MDF refers to the stellar particles in the simulation which are counter-rotating with 
v 0 < 0 at a projected radius R > 15 kpc. 
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Figure F.5: Kinematically selected halo MDFs for the M to t = 5xl0 12 M 0 semi- 
cosmological simulations in Renda et al. (2005b). The 68% Confidence Level range 
and the number of stellar particles each MDF relates to are also shown. Each 
MDF refers to the stellar particles in the simulation which are counter-rotating 
with v 0 < 0 at a projected radius R > 15 kpc. 
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Figure F.6: (continued) Kinematically selected halo MDFs for the M tot = 5xl0 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF relates to are also shown. Each 
MDF refers to the stellar particles in the simulation which are counter-rotating with 
v 0 < 0 at a projected radius R > 15 kpc. 
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Figure F.7: (continued) Kinematically selected halo MDFs for the M tot = 5xl0 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF relates to are also shown. Each 
MDF refers to the stellar particles in the simulation which are counter-rotating with 
v e < 0 at a projected radius R > 15 kpc. 
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Figure F.8: (continued) Kinematically selected halo MDFs for the M tot = 5xl0 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF relates to are also shown. Each 
MDF refers to the stellar particles in the simulation which are counter-rotating with 
v 0 < 0 at a projected radius R > 15 kpc. 
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Figure F.9: (continued) Kinematically selected halo MDFs for the M tot = 5xl0 12 M 0 
semi-cosmological simulations in Renda et al. (2005b). The 68% Confidence Level 
range and the number of stellar particles each MDF relates to are also shown. Each 
MDF refers to the stellar particles in the simulation which are counter-rotating with 
v e < 0 at a projected radius R > 15 kpc. 
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APPENDIX G. STELLAR HALO (V - I) COLOUR DISTRIBUTIONS 



1 1.5 2 2.5 3 3.5 1 1.5 2 2.5 3 3.5 

(V - I) 

-3.5 < I < -3.0 (v - i) 

-3.5 < I < -3.0 


Figure G.l: Halo synthetic CMD (V - I) colour distribution at -3.5 < I < -3.0 for 
the M tot = 5xlO n M 0 semi-cosmological simulations in Renda et al. (2005b). Each 
panel also shows the number of generated stars at -3.5 < I < -3.0 in each synthetic 
CMD. 
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Figure G.2: (continued) Halo synthetic CMD (V - I) colour distribution at - 
3.5 < I < -3.0 for the M tot = 5xlO n M 0 semi-cosmological simulations in 
Renda et al. (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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Figure G.3: (continued) Halo synthetic CMD (V - I) colour distribution at - 
3.5 < I < -3.0 for the M to t = 5xlO n M e semi-cosmological simulations in 
Renda et al. (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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Figure G.4: (continued) Halo synthetic CMD (V - I) colour distribution at - 
3.5 < I < -3.0 for the M tot = 5xlO n M 0 semi-cosmological simulations in 
Renda et al. (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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Figure G.5: Halo synthetic CMD (V - I) colour distribution at -3.5 < I < -3.0 for 
the M tot = lxlO 12 M 0 semi-cosmological simulations in Renda et al. (2005b). Each 
panel also shows the number of generated stars at -3.5 < I < -3.0 in each synthetic 
CMD. 
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Figure G.6: (continued) Halo synthetic CMD (V - I) colour distribution at 
-3.5 < I < -3.0 for the M tot = 10 12 M 0 semi-cosmological simulations in 
Renda et ah (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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Figure G.7: (continued) Halo synthetic CMD (V - I) colour distribution at 
-3.5 < I < -3.0 for the M to t = 10 12 M 0 semi-cosmological simulations in 
Renda et al. (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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Figure G.8: (continued) Halo synthetic CMD (V - I) colour distribution at 
-3.5 < I < -3.0 for the M tot = 10 12 M 0 semi-cosmological simulations in 
Renda et ah (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 




196 


APPENDIX G. STELLAR HALO (V - I) COLOUR DISTRIBUTIONS 


3 

2 

1 

0 

3 

I 2 

1 

0 

3 

i 2 

i 

0 

Figure G.9: (continued) Halo synthetic CMD (V - I) colour distribution at 
-3.5 < I < -3.0 for the M to t = 10 12 M 0 semi-cosmological simulations in 
Renda et al. (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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Figure G.10: Halo synthetic CMD (V - I) colour distribution at -3.5 < I < -3.0 for 
the Mt 0 t = 5xl0 12 M 0 semi-cosmological simulations in Renda et al. (2005b). Each 
panel also shows the number of generated stars at -3.5 < I < -3.0 in each synthetic 
CMD. 
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Figure G.ll: (continued) Halo synthetic CMD (V - I) colour distribution at - 
3.5 < I < -3.0 for the M to t = 5xl0 12 M 0 semi-cosmological simulations in 
Renda et al. (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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Figure G.12: (continued) Halo synthetic CMD (V - I) colour distribution at - 
3.5 < I < -3.0 for the M to t = 5xl0 12 M e semi-cosmological simulations in 
Renda et al. (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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Figure G.13: (continued) Halo synthetic CMD (V - I) colour distribution at - 
3.5 < I < -3.0 for the M to t = 5xl0 12 M 0 semi-cosmological simulations in 
Renda et al. (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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Figure G.14: (continued) Halo synthetic CMD (V - I) colour distribution at - 
3.5 < I < -3.0 for the M to t = 5xl0 12 M e semi-cosmological simulations in 
Renda et al. (2005b). Each panel also shows the number of generated stars at 
-3.5 < I < -3.0 in each synthetic CMD. 
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APPENDIX H. HALO MDF VIA METALLICITY-COLOUR 
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Figure H.l: Halo MDF via Metallicity-Colour relationship by the same pipeline as 
in Mouhcine et al. (2005) for the M to t = 5xlO n M 0 semi-cosmological simulations 
in Renda et al. (2005b). The 68% Confidence Level range and the number of stellar 
particles each MDF refers to are also shown. 
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Figure H.2: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M to t = 5xlO n M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure H.3: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M tot = 5xlO n M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure H.4: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M to t = 5xlO n M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure H.5: Halo MDF via Metallicity-Colour relationship by the same pipeline as 
in Mouhcine et al. (2005) for the M tot = 10 12 M 0 semi-cosmological simulations in 
Renda et al. (2005b). The 68% Confidence Level range and the number of stellar 
particles each MDF refers to are also shown. 
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Figure H.6: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M tot = 10 12 M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 











210 


APPENDIX H. HALO MDF VIA METALLICITY-COLOUR 



-4 -3 -2 -1 0 -4 -3 -2 -1 0 

[Fe/H] [Fe/H] 


Figure H.7: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M tot = 10 12 M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure H.8: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M tot = 10 12 M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure H.9: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M tot = 10 12 M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure H.10: Halo MDF via Metallicity-Colour relationship by the same pipeline as 
in Mouhcine et al. (2005) for the M to t = 5xl0 12 M 0 semi-cosmological simulations 
in Renda et al. (2005b). The 68% Confidence Level range and the number of stellar 
particles each MDF refers to are also shown. 
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Figure H.ll: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M tot = 5xl0 12 M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure H.12: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M to t = 5xl0 12 M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure H.13: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M tot = 5xl0 12 M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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Figure H.14: (continued) Halo MDF via Metallicity-Colour relationship by the same 
pipeline as in Mouhcine et al. (2005) for the M to t = 5xl0 12 M 0 semi-cosmological 
simulations in Renda et al. (2005b). The 68% Confidence Level range and the 
number of stellar particles each MDF refers to are also shown. 
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